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Abstract 

The present review includes the description of theoretical methods for the inves- 
tigations of the spectra of hydrogen-like systems. Various versions of the quasipo- 
tential approach and the method of the effective Dirac equation are considered. 
The new methods, which have been developed in the eighties, are described. These 
are the method for the investigation of the spectra by means of the quasipotential 
equation with the relativistic reduced mass and the method for a selection of the 
logarithmic corrections by means of the renormalization group equation. The spe- 
cial attention is given to the construction of a perturbation theory and the selection 
of graphs, whereof the contributions of different orders of a, the fine structure con- 
stant, to the energy of the fine and hyperfine splitting in a positronium, a muonium 
and a hydrogen atom could be calculated. 

In the second part of this article the comparison of the experimental results 
and the theoretical results concerning the wide range of topics is produced. They 
are the fine and hyperfine splitting in the hydrogenic systems, the Lamb shift and 
the anomalous magnetic moments of an electron and a muon. Also, the problem 
of the precision determination of a numerical value of the fine structure constant, 
connected with the above topics, is discussed. 
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1 Introduction 



In the last years, investigations of energy spectra of hydrogen-like systems, as a positronium, a 
muonium, a hydrogen atom e. t. c., are in progress. These two-particle systems can be used as a test of 
quantum electrodynamics (QED). The former review [El has been devoted to the consideration of this 
problem. The modern status of investigations in the above direction is analyzed in the present paper. 

It is well-known, the Dirac equation gives the opportunity to take into account the spin-orbit in- 
teraction in hydrogenic systems and, correspondingly, to predict the fine structure of energy levels with 
the accuracy 0(a ). However, the changes in the fine structure levels, the corrections to the hyperfine 
splitting (HFS), the Lamb shift (splitting of the 2S , 1 / 2 an d 2Pi/ 2 levels which coincide in the Dirac theory) 
can be described on a basis of quantum field methods only. 

In nonrelativistic quantum mechanics, the two-body problem is reduced to simpler ones, motion of 
the center of mass (c. m.) and motion of a particle with a reduced mass. In the relativistic case, QED, 
the separation of c. m. motion is impossible. Also, definition of such a notion as a potential is impossible 
in a usual way. The remarkable feature of quantum field methods for the description of bound states is 
the use of the formalism based on the two-particle Green functions. Then, the spectrum is found as the 
positions of its poles. 

The equation for the two-fermion Green function could be written in the Schwinger form J2|: 

{( 7 7r-M) 1 ( 7 7r-M) 2 -7 12 }G = 7, (1.1) 

where 7r a = p a — eA a with p l a being 4-momenta of i particle, A l a is an external field influencing i particle; 
e is an electron charge; Mi is a mass operator for i particle, Ii 2 is an interaction kernel between 1 and 2 
particles; I is a unit operator. Finally, 

G(x 1 ,x 2 ;x 3 ,x,) = <Q|5|0> (L2) 

is the total two-particle Green function in the interaction representation, ip(xi) are the field operators of 
constituent pa rtic les. 

Equation ( p~l| ) could be re-written in the Bethe-Salpeter (BS) form 

G(xi, x 2 ; x 3 , X4) = G (xi,x 2 ; x 3 , x 4 ) + G (xi,x 2 ; x' 1 ,x' 2 )K B s(x' 1 ,x' 2 ; x' 3 , x'^Gix^x'^ x 3 , £4), (1-3) 

G a (xi,x 2 ;x 3 ,x i ) = iG a (xi,x 3 )G b (x 2 ,X4,), (1.4) 

where G a ^ b are the Green functions of free fermions, K B s is the kernel for the BS equation, which is 
connected with the interaction operator of particles I12 which is a sum of the two-particle irreducible 
Feynman diagrams. The state of the two-particle system is defined by the two-time wave function tj), 



which is a solution of the homogeneous equation corresponding to (1.3) 



(Go" 1 -K BS )*v{x 1 ,x 2 ) = 0, (1.5) 

^ v {x u x 2 ) =< I T{Mxi)Mx2)} \ V,v>. (1.6) 

The ket- vector | V, v > characterizes the bound system as a whole with the 4 - momentum V and a 
set of additional quantum numbers v. 

Using a translation invariance and choosing the c. m. system 7^= (E, 0) one can obtain the wave 
function. It corresponds to the state of the definite value of energy E 

■*v(x 1 ,X2)=e- lEXa dp E (x), (1.7) 

Xq is the time component of the cm., x are the relative coordinates. 

The bound state problem can be solved in the relativistic quantum theory only approximately, by 
means of the perturbation methods. The primary approximation is usually chosen so that it corresponds 



1 Integration is implied in (1.3) and further on the repeated variables. 
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to the instantaneous Coulomb interaction. Then, an energy spectrum consists of the Coulomb levels 
found from the wave equations. The corrections to the energy spectra are obtained from the higher 
orders of perturbation theory Gj 

AE = -i$ Ko (x){K + KG c K + e. t. c. )4>k c (x'), (1.8) 

where K = Kbs — Kc , Kc is the Coulomb part of the kernel of the BS equation, Gc is a solution 



of Eq. Q1.3Q with the kernel Kc, 4>k c ( x ) is a solution of Eq. (1.5) with the kernels Kc- However, the 
state 4>k c i x ) is time - dependent; the relation of the function 4>k c with a solution of the Schrodinger or 
Dirac equations with the Coulomb kernel is sufficiently complicated. The normalization and formulation 
of the boundary conditions are not clear for the wave function depending on the relative time. All the 
above-said influences the calculation accuracy. 

The formalism based on the three-dimensional equations has been produced in relativistic bound 
state theory even before appearance of the covariant formalism in quantum field theory Q — H. In this 
connection, an important meaning did have the development of the quasipotential approach |g, ^] and 
Gross method [jl0|. In these approaches quantum field equations are deprived of shortcomings of the 
BS equation and they are formally close to the nonrelativistic Schrodinger equation with a potential V . 
They can be considered as direct generalizations of the potential two-body theory to the relativistic case. 
The principal idea of these three-dimensional methods consists in a choice of a " primary" two-particle 
propagator which describes well physical meaning of problem. In the quasipotential approach, this choice 
is made by the conversion to the two-time Green function in which the parameter of a relative time t a — tb 
is equated to zero. In the Gross approach, the "primary" two-particle propagator is chosen as a projector 
onto the mass shell of a heavy particle and an electron propagator. 



2 Quasipotential approach in quantum field theory 

The quantum field equations for a system of two particles, mentioned in the Introduction, are 
reduced in the three - dimensional approaches to equations like the Schrodinger ones with the quasipo- 
tential defined by the two-time Green function. In spite of the lack of a clear relativistic covariance, the 
quasipotential method keeps all information about properties of the scattering amplitude which could 
be received starting from the general principles of quantum field theory. Therefore, one can investigate 
both analytical properties of the scattering amplitude, its asymptotic behavior and some regularities of 
a potential scattering, e. g., at high energies [|ll| . The renormalization procedure of the quasipotential 
equation is reduced to a charge and mass renormalization like the usual S- matrix theory Jl^| . 

The quasipotential method is also very efficient for determination of the relativistic and radiative 
corrections to energy spectra of hydrogen-like atoms. In some cases, it is convenient to define the two- 
particle off-shell scattering amplitude instead of the Green function (ll. 



G — Go + GoTGq, (2-1) 
which is connected with the kernel of the BS equation by the expressions: 

T = K BS + K BS G T (2.2) 

or 

T = K BS + K BS GK B s- (2.3) 



On the mass shell (po — qo — 0, \Jp 2 + rnf + \/p 2 + m\ =j_ \/cf 2 + m\ + y/q 2 + m| = E) the amplitude 
T is equal to the physical scattering amplitude (see Fig.lp. 

The quasipotential approach is universal and symmetrical in describing both the particles. Due to 
this, it is used for consideration of any system of arbitrary mass particles. However, the way of equating 



2 The physical sense of the quantities po,qo, p, q is clear from Fig. 1; mi and ttt-2 are masses of the 
constituent particles. 
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the times does not permit one to get the Dirac equation for one particle when mass of another tends 
to infinity. When we consider hydrogenic atoms with large nucleus charge and large mass, it would be 
convenient to use (as the initial approximation) the exact solution of the Dirac equation with the Coulomb 
potential. This procedure has been proposed in J[o| |. 

In the paper ]T3], the three-dimensional self-according formalism has been built, which leads to the 
effective (modified) Dirac equation, following to the ideas of the quasipotential approach. In contrast 
with the usual procedure of equating the times in the operators ip a (%i) and ipb(x2) of Eq. (1.6), the 
procedure of tending x\ — > 00 to infinity a has been carried out in the above paper. After this, the 
following expression has been obtained for the WF in the momentum representation: 



<P2,&2 I V-a(O) \ V,V>, 



(2.4) 



in which it is easy to go over to the variables of the total and relative momenta. A similar procedure, 
carried out for the free Green function Go ( see Eqs. (2.5), (3.4), (3.5) of the article (l3|), leads us 
directly to the equation: 



(j]\V +P~ mi)^) v (p) = 



1 



(27T) 



dq- V(p,q)^ v (q), 



(2.5) 



which as 7112 — * 00 goes over in the Dirac equation for the first particle in an external field. This property 
is a characteristic feature of the method of the effective Dirac equation (EDE). 

Although the Logunov-Tavkhelidze quasipotential approach and the Gross approach (the EDE method) 
have different physical foundation, the formulae used to calculate the energy level corrections in compos- 
ite systems are very similarcl. 



Quasipotential approach 



EDE Method. 



The initial 
for a two 



propagation function 
fermion system: 



G Q {p,q;E) = 

OO OO 

= (2^F / / d P° d( loGo(p,q;po,qo;E); 

Go{p,q;p ,qo;E) = 

= 1 (2k) A ^—^—8^ (p - q). 

v ' pi-mi p2 —m'2 v ' 



g(p,q;E)=A 2 S(p,q;E) = 

= 2 m e(p Q 2 )S(pl- m l)(2n) A m±^S^(p-q). 



The relative time parameter in this 
approach is equated to zero. In the 
momentum representation integration over 
the relative energies corresponds to this 
procedure. 



A2 denotes the projector onto the mass 
shell of the second particle. The momenta of the 
components are usually written as p\ = V — P2 ■ 
In the case of the c.m.s., namely V — (E,0), 
one has p\ = E — y/p 2 + mf. 



Gq = F = {2nfS(p- q)(E - e lp - ^p)- 1 
is a free particle Green function projected 
onto the states of positive energies. 



3 When we consider the Green function (1.2), it is also necessary to carry out the operation x\ — ► 00. 
4 The definitions of momenta correspond to Fig.l. 
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The total Green function. 



The total two-time Green function 
is connected with the off-shell scattering 
amplitude: 



G+ = G+ 



The effective Dirac equations is built 
by means of the new two-particle Green 
function G: 

G = g + g-T-g, 



which can be re-written using the quasipo- 
tcntial 



t> = (G+)" 1 



(G+)- 1 



(Go+r 1 

= T 



GqTGq 



1 + G+rj , 
in such a form: 

G+ = G+ + G+ -V -G+. 



which satisfies the equation: 

G = g + g- Vede ■ G. 

By means of this Green function the three- 
dimensional Green function G can be 
defined: 

G(p,q;E) = 

= i(2nfS(pl - m 2 2 )G(p, q; E)S(q 2 2 - m|), 

which satisfies the analogous equation. 

The equation kernel 



is the quasipotential expressed by means of 
the amplitude 



r = (G+)- 1 [GoTGo 
as follows: 

t = V + VG+t- 



is introduced by the equation: 
T = Vede + Vede ■ 9 • T. 

The wave function equation. 



The single-time wave function obeys 
to the corresponding homogeneous 
equation: 



The homogeneous equation for the 
wave function in a symbolic form is 



(G+)- 1 * = T>*. 



S- V = A 2 • Vede 
The expanded form of which are 
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(pi - m-i)<j)(j); E) = 

= (p 2 + m 2 ) J {2 irp2E q WEDEiP, f, E)(/>(q; E). 

Since the WF on the muon index satisfies 
the Dirac equation, the above mentioned 
equation is re- written for the WF 'J, having 
the electron spin factor only: 

CP! - mi)f(p; E) = J j^iV EDE {p, q; E) 



The energy spectrum. 

After solving the obtained equations by means of perturbation theory the formulae for the corrections 
to energy levels of a hydrogen-like atom are |Q 

E n = E°+ < n | AV^ + + ... \ n > x E n = E° n + < n \ iSV ED E \ n > x 

x(l+<n\^^\n>)+ U +<n \ l dSVm>B\ n> \ + 



8E I J I 8E 

+ <n\ E m# „ AVWl^^gAV^ | n > +... +<n\ iSV EDE G nQ i6V EDE \n>x 

x(l+<n\ i a % E E DE | n >) + ... 
AV^ = — vc, Vc is the Coulomb potential. Vede = Vq + SVede, Vq is the "primary" potential. 

(The derivative on E is taken at the point E^.) 




Fig. 1. The parametrization of the two-particle off-shell scattering amplitude T in c. m. s. 

i 7 ?! 2W 1 ' r ' 2 2m )• 

It should be noted that there exists another way of constructing the quasipotential by using the 
physical on-mass-shell scattering amplitude but not on energy shell. However, it is extremely important 
to take into account the binding effects and relativistic interaction effects on precision calculations of the 
eigenvalues of hydrogen- like atoms, e. g., on calculations of the corrections to the Fermi energy of the 
HFS of the ground state with the accuracy higher than a 5 . The efficiency of an analysis of these effects 
is essentially greater in the framework of the first method. The cause is that the on-shell scattering 
amplitude can possess singularities in the infrared region in higher orders of perturbation theory. If we 
are concerned with the diagrams up to the fourth order of charge, the infrared divergencies cancel one 
another at the end. In analysing higher orders in a, the singularities of the quasipotential lead to the 
logarithmic contributions with respect to the fine structure constant. The possibility of including these 
effects on the basis of the on-shell scattering amplitude is restricted. Moreover, the existence of poles 
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in the virtual particle propagators is a complicated factor in integrating over the relative momenta (see 
Section III and Q). 

When the quasipotential is built by the first method on the basis of the two-particle Green function, 
the total energy appears in the quasipotential directly. As a result, the normalization condition and the 
condition of orthogonality of the eigenfunction Jil| and at the same time perturbation theory acquire a 
specific peculiarities. 

Let the equation of the eigenvalue problem has the form: 

[F-^E) - V(p,q;E)] V E (q) = 0, (2.6) 

wherei 

F(E) = G+ = (2ir) 3 6^(p - q)(E - e lp - e 2p )-\ (2.7) 

and the quasipotential be defined by Eq. on page No. 4. Then, we have the following expression for the 
eigenfunction ^> n corresponding to the state with the eigenvalue E n : 

G + (E)W(E, E n )* n = — , (2.8) 

with 

W(E, E.) - ^<*> - M , m 

hi — £j n £j — £j n 

Since near the singularity (E ~ E n ) one has 

6 + (E)~-^L., (2.10) 

the orthonormalization condition has the form: 

V* m W(E m ,E n )y n = 6 mn , (2.11) 

and 

dV(E) 

4-™*„ = S mn + ^tn^T 1 E=E n {E m = E n ). (2.12) 

oE 

Since the energy levels of a composite system are the poles of the exact scattering amplitude and 
every separate term of a charge power expansion of the amplitude does not have these poles, there must 
exist an infinity series of the diagrams for bound states the contributions of which to the energy of a 
bound system have the same order on a. The diagrams, which this series consists of, are the reducible 
two-particle Feynman diagrams. However a similar binding effect takes place also for the irreducible 
diagrams, e. g., when calculating the Lamb shift. 

The diagrams of successive Coulomb photon exchanges are shown in ]3|, [lTj to contribute the same 
order in a in the infrared region when calculating the HFS of hydrogen-like systems. So, when constructing 
the quasipotential, the corresponding modification procedure is necessary because we must sum the 
infinite series of diagrams selectively. Let us introduce the Coulomb Green function Gc satisfying the 
equation 

(G Q 1 -K c )Gc=I, (2.13) 

where 

e(Ze)T a 
(p ~ q) 2 

is the Coulomb kernel (L = 7io72o)- The expression ( p,13| ) can be rewritten in the form: 



K C (p, q) = ~ ;\ = v c To (2.14) 



Gc = Go + GoKcGc = Go + GcK c Go, (2.15) 



5 In view of the absence of the inverse Green function to the Green function of the system of free two 
fcrmions, it is necessary to carry out the procedure of projecting Go onto the states of positive energies 
in the quasipotential approach. 
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and the total Green function can be represented similarly to (p.l|- |2.3| ): 

G = Gc + Gc'TGc, 

f = K + KG C f, 
f = K + KGK, 



(2.16) 

(2.17) 
(2.18) 



with K = Kbs ~ Kc- 

The wave function specifying the state of the definite energy E of the two-fermion system must satisfy 
the quasipotential equation 



(G+(P, T, E))- 1 - V(p, q; E) * E (q) = 0, 



(2.19) 



where once again an integration is implied on the repeated momenta. The corresponding quasipotential 
has the form: 



V(p,q:E) = G+(p,q;E) 



G + (p,q;E) 



{G + C Y 



GcTGc 



(2.20) 



and G + is the two-time Green function projected onto the positive energy states (see above). 

The technique developed in |Tg|] , which uses the Coulomb Green function formalism, gives us the 
opportunity to take into account the successive multiplied exchange by the Coulomb photons and to 
write formulae for the energy level shift AE with respect to the Coulomb one. 

From Eq. ( [2.16 ) it is clear that the total two-particle Green function (projected onto the states of 
positive energies) has the inverse one 



(G + )- 1 = (G+)- 1 -(G+)- 1 f+(G+)- 1 , 



i c — 



GcTGc 



and the quasipotential can be presented in the form: 

V(p,q;E) = (G+)- 1 f+(G+)- 1 =f c - fcG+fc + 

f c = (G+)" 1 f+(G+)- 1 
After the mentioned projecting of the Coulomb Green function, one can find 

G+=F + F{v c G c ) + =F+ {G c v c ) + F. 

Here, the expression has been used: 

Go = (A++F-A— F')T , 



where F is defined by Eq. (2.7) and 

F' = (2n) 3 8(p- q)(E + ei p + e 2p y l . 

A ++ = A^ (p)A2 (— p), A = Aj~(p)A^"(— p) , A ± (p) are projecting o perat ors. 
The following formulae are used for some transformation of Eq. (|2.24 ): 



[QGcY 
[GcQ]' 



[Qsr ] + G+, 
G+ [ST Q] + , 



(2.21) 

(2.22) 
(2.23) 

(2.24) 
(2.25) 
(2.26) 



(2.27) 
(2.28) 



where S = (1 + A F'vc) , £ = (1 + ucA F') 1 , and Q has an arbitrary matrix structure. As a 
result, we get the closed equation for the function Gjt: 



Gj = F + FK-zGq = F + G^K^F 



(2.29) 
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with the kernel As = w|u2^ u i M 2, S = t>c£ — S'uc- In this case, the inverse Green function (Gq) 1 is 



(Gj) = .F — K-s (2.30) 
and the kernel K s includes the projections of a Coulomb interaction onto the negative energy states 



= A" 



c 



(2.31) 



Let us note that F' , in contrast with F, has no singularity when p, q — > 0, E — > mi + ma. In the 
approximation of the upper components (iti = (™'),Wj is a normalized Pauli spinor.) the kernel Kq 
coincides with the Coulomb potential and the term, including the projecting operator A , is equal to 
zero. Thus, the principal part of the kernel As is equal to Kq. The function (Gjfc) -1 coincides with the 
inverse Green function of the non-relativistic Schrodinger equation with the Coulomb potential. 



In many cases, the quasipotential is conveniently expressed by the amplitude Tq 
can be simply carried out by the helpful formulae ( [2.27| .2.28). Since 

Tc = T o + (f K c G c ) + + (G c K c f ) + + {G c K c f K c G c ) + , 



GqTGq 



we have after transformations (2.27, 2.28) : 

f+ = f + + G+(S'« c fo)+ + (f oUc r i]r )+G+ + G+(S'« c f A' c Sro) + G+. 



The expression (2.33) is re- written in a more convenient form: 

Tq = (GquIu^ ■ I + I ■ U*U2)To(TqUiU2 ■ I + Tq ■ IYiUiU2Gq), 

or, with taking into account 

\+ + \-=I, X- =A+(p)A 2 -(-p) + Ar(p)A+(-p) + A— , 

in the form: 

f + = Gj If^t+f- 1 + F- 1 (f r a \-±r ) + + (ibx-f^+F- 1 + (£A-T r A-£ro) + 



G 



This 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



Therefore, the amplitude tq can be presented as the sum tq — f, o + p, where fo is the scattering 
amplitude including the T- matrix. The rest of the terms corresponds to interactions of higher order with 
subtraction of iterations. 

The graphical interpretation of the amplitude fo is clear. It is a series of all irreducible diagrams with 
the exception of the one-Coulomb exchange diagram and the reducible diagram with Coulomb exchanges 
in intermediate states only. The corresponding quasipotential can be divided in two parts: 



V(p,q;E) = V fo +V p 



higher orders of perturbation expansion, 



where 



with 



V„ 



"T — fo - foG^fo, 
P - rGqP - T G^p - pG^T , 



p = F^ARF' 1 , 
AR = F(t\-f ) + + (Tbr A-Sr ) + A + A(£A-T r A-£ro) + A. 



After substitution of the amplitude T from (2.18) (T ~ K + KG C K) into (2.38) we have: 



Vf 



F- 1 



K+ + (KG C K)+ - K+F- 1 G+F- 1 K+ + e. t. c. 



F~ 



(2.37) 

(2.38) 
(2.39) 

(2.40) 
(2.41) 

(2.42) 



s 



in which the notation is used: 



(KG C K)+ = 



GqKGq 



GqKGcKGq 



(2.43) 
(2.44) 



For calculations of the energy level shift up to the accuracy a e log a, it is sufficient to take into account 

{KG C K)+ = (KG k)+ + (K Q K c ko) + + {K T Q v c G c v c T Q K ) + . (2.45) 

Using the representation of a unit operator through projecting operators we can find out 

(K Q roIv c Gcvcr lko) + = K+{v c G c v c ) + k+ + K+{v c G c v c To\~K ) + + 

+ {KvT \-v c G c vc) + k+ + {K r \-v c Gcv c r \-Ko) + . (2.46) 

The first term of the expression derived contributes mostly in calculating the energy spectra. 
Using the definition of the operator E, the following transformation can be made 

(v c G c v c ) + = [«oS(J + K—F'vc)G c vc\ + = (2.47) 

= {±\+G c vc) + + {t\-G c v c ) + + (EA— F'v c Gcvc) + . (2.48) 
Starting from the definition of the kernel Kj; ( [2.31 ) , we get 

(t\+G c vc) + = Kv(G c v c ) + = KzG+Kv (2.49) 

Since the equality \~Gc = — A~~F'(ro + vcGc) holds, the sum of two last terms of the expression 
( 2.48 ) can be transformed in 



Therefore, 



(£A-G c «c) + + (EA~ F'v c G c vc) + = -(EA~ F'I> C ) + = 



(v c G c v c ) + = K^G+Kv +Kz-K+. 



(2.50) 



(2.51) 



The iteration term in the quasipotential expression ( 2.42| ) is written in the following way using ( 2.29| ): 

K+F- 1 G+F- 1 K+ = K+iF- 1 +K^+ K^G+K^)K+. (2.52) 

Substituting the last two formulae in the quasipotential expression we can be convinced that the pro- 
jection of the block vcGcVc onto the positive- energy states drops out from the quasipotential expression 
fl2.42p . In this case 

Vf = F- 1 {k+ + (KG a K)+ - K+F- 1 ^ + (k Q K c k )+- 



K+K+K+ + k+(v c Gcv c T \-k ) + + {k T \-v c G c vc) + k+ + 



{k ro\-vcGcvcr \- k ) + ] f~\ 



(2.53) 



The information about the corrections to the Coulomb levels c an be obtained by constructing per- 
turbation theory on the basis of the above presented ideas and Eq. ( 2.19 ) in the form: 



(f-\E c ) + AE-K+- U(£)) V E = 0, 



where 



V = V + Kv-K? 



(2.54) 
(2.55) 
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AE = E — Ec is the correction to the ground state of the non-relativistic Schrodinger equation with the 
Coulomb potential. 

If E' is an eigenvalue of the wave function of the equation with the kernel Kq , 

(G+(E'))-^e' = (F- l (E') - K+) * E , = 0, (2.56) 



the initial equation ( P.54|) gets the following form: 

((G+OEOr 1 + AE' - V(E)) * £ = 0, (2.57) 

with AE' = E - E' and AE = AE' + AE C = AE' + E' - E C - Let us assume that the eigenfunction 
^ E' and the eigenvalues are known. In accordance with the methods of perturbation theory in the first 
order we have e = ^ E' + 

((G+(£')r 1 1'/ + AE'^e' - V(E')^e') = 0. (2.58) 

Having multiplied both of parts of this equality by ty E , on the left side and using the normalization 
condition we come to: 

AE'j = < V E > | V{E') | ^ E > > (2.59) 

= (g+(E')-^§-)v(E')^e>, E^E'. (2.60) 



E - E' 

Analogously, for the correction to the energy level in the second order we have 

AE'jj =< V E > | V(E') (l + G' c (E')V(E'j) \ * E < >, (2.61) 

where 

G'UE')=G+(E')-^§-. (2.62) 

The term including the energy derivative disappears from the final result because ^> e' is the eigenfunctions 
of the eqution with the kernel K. % wh ich does not depend on the energy. 

To determine \1/ e 1 from Eq. ( 2.56| ) let us introduce the auxiliary function $ e 1 defined by the equality 

2fi(E + ei q + e 2q ) 

(2E'e lq + E' 2 +m 2 - m 2 )(2E'e 2q + E> 2 + m| - mf) 

= a vi 21 vi i i \ E ' W> [2.b6) 

8/j.E' Z (E' + ei q + e 2q ) 

which satisfies the non-relativistic Schrodingcr-likc equation 

(g^ 1 - Ae - 8K C )$E' = 0. (2.64) 

P 2 

9 C = w c - y~ - v c , Wc=E c -mi-m 2 , (2.65) 

((mi + m 2 ) 2 - E' 2 ) (E' 2 - (mi - m 2 ) 2 ) 
Ae = W c + ^ - — - 11 =s -AE C (2.66) 



Here 



and 



8K C = K+(p, - Vc ^ K+-— - v c . (2.67) 



Xe i q + mE){t 2q + m E) _ K+ W c 

2»(E + e lq + e 2q ) V ° c E> - e lq - e 2q 
The Coulomb WF describing the IS* state is: 

^= ^r/w) 2 iMr=o)i ' (2 - 68) 
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|^( r = )| 2 =M! (2.69) 

It satisfies the equation 

9c Vc - 0. (2.70) 

The correction to the Coulomb energy level AEc and the function $£/ are found (as above AE') by 
means of a quantum-mechanical perturbation theory 

AE C =< <i>c I SK c (l + g c SK c ) \ <f> c >, (2.71) 

9c =9c- e C -\ c ' when ^^- ( 2J2 ) 



In the second order of perturbation theory the eigenfunction of Eq. (2.56) has the following form: 
= (e lp + mi )(e 2p + m 2 )-f r + . ^ . i 

4mi77T,2 L J 

(Integration is implied over 3- vectors k and <f.) 

The final expression for the total shift of energy levels with respect to the Coulomb level is the sum of 
the correction AEc and the corrections from the quasipotentials SKs — K-£ — Kq and V(E). In [n| 
the following expressions for these corrections are presented: 

AE = AE KK + AE KV + AE VV , (2.74) 

AE KK = < 0c I SK a (l + g'JKc) \(f>c> + 

+ < V B , | SK S (1 + G' C (E')SK S ) | * B > >, (2-75) 

AE KV = < i$E> | f (£')G C (£')<^£ +SK S G' C (E')V(E') \ >, (2.76) 

AEvv = <^E'\ V(E')(1 + G' C (E')V(E')) \V E '>- (2.77) 

The problem of construction of the kernel for the quasipotential equation was under consideration in 
p!^-p0|. The kernel K can be expanded in the series of perturbation theory: 

K = K T + K ( - 2 \ (2.78) 

where the index T denotes the transverse photon in the Coulomb gauge; is built from the diagrams 
of the second order in the fine structure constant. 

Let us point out that when the interaction kernel for the two-fermion system including a particle and 
an antiparticle is constructed, it is necessary to take into account the annihilation interaction channel 

To investigate the energy spectra of a two-particle relativistic bound system, the method has been 
proposed which is based on the use of the local quasipotential equation with the relativistic reduced mass 
in the c. m. s. 124]. By means of the "rationalization" of the Logunov - Tavkhelidze equation |^5| 

the following equation has been obtained: 

( b^E)__JP_ E I d? v(fitE)i6E( fi (2J9) 

2fi R 2fi R J (2tt) j 

where 

I{E,p) = {E + £lp + e2p){E2 (eip £2p)2) . (2.80) 
&EiE 2 E 

In accordance with the definitions of Fig.l: 

E 2 — m| + m? 

Ei = mE = 2J K (2.81) 

E 2 — m? + 77X9 

E2 = mE = 2E ' (2 - 82) 
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Ei + E 2 = E is the energy of a bound state; 



b 2 (E) = ^L[E 2 - (mi + m 2 ) 2 ] [£ 2 - (mi - m 2 ) 2 ] (2.83) 

is the squared relative momentum on the energy shell. Finally, 

E X E 2 E X E 2 E 4 - (m 2 - m 2 ) 2 

— - e7Te- 2 4iP (2 ' 84) 

is the relativistic reduced mass defined in accordance with the relativistic expression of the coordinate 
of the center of mass. Let us mention other definitions of the relativistic reduced mass like jiR = 
\ E2 -^- m ^ 2 ] ; Ref H ^/ = m^ii^ Ref H The quant it y m > = ^fmTm^, the effective mass, which 
has been defined in Ref. |27[ ], is closely connected to the notion of the relativistic reduced mass. This 
definition gave the authors of |^7j the opportunity to reduce the relativistic two-body problem to the 
case of the particle motion with the mass m! in the quasipotential field. In the non-relativistic limit 
E\.2 — » tii, 2 the relativistic reduced mass /ir becomes the well-known reduced mass fj, = ™ i 1 J^ 2 ■ 

As has been pointed out above, the quasipotential can be considered from the Lippmann-Schwinger- 
like equation with the scattering amplitude with the relative energies of particles pa = qo = being equal 
to zero, pi = Ex, p 2 = E 2 . On the energy shell E = e ip + e 2p = ei q + e 2q , e lp = p 2 + m 2 , e lq = 
\Jq 2 + rnf we have I(E,p) = 1, and in this case the equation (2.79) can be solved exactly for the 
Coulomb interaction. 

In the case of an interaction of two spinor particles with masses mi , to 2 and charges (-e) and (Ze) the 
main contribution to the binding energy of particles is shown in p2[ | to arise from the modified Coulomb 
potential 

^^) = --^(l + ^|). (2.85) 
Quantization of energy levels is defined by the equality analogous to the equation obtained in p8[ : 

b 2 E 2 Z 2 a 2 



{b 2 + ExE 2 ) 2 



(2.86) 



(n = 1, 2 ... is the principal quantum number). The above equality leads to the variant of the relativistic 
Balmer formula p9| : 

(Za) 2 ^- 1/2 



E z n = m{ + mi + 2miTO 2 (^1 + j . (2.87) 

The above formula can be re- written as an expansion of the binding energy B in powers a 2 : 

n 7? nZ 2 a 2 fiZ 4 a 4 [i 2 m, x m 2 

B = E-mx-m 2 = -- — — + - — — (3 ), m= ■ (2-88) 

A n z on 1 m\m 2 m\ + m 2 

The relativistic Balmer formula takes into account the recoil effects but does not describe the fine and 
hyperfine structure corresponding to the spin-orbit and spin-spin interactions. These corrections have 
been considered in |24| using the formalism of the local quasipotential equation with the relativistic 
reduced mass. In the first order of perturbation theory we have 

AEi =< *' c | AV y + V 2y + e. t. c. | >, (2.89) 
where AVy = V 7 — V™ is determined by the difference of the one-photon exchange quasipotential and 



the modified Coulomb potential (2.85). In turn 



T/ T+f*- n M f dk V~f(P^;E)V 7 (k,q;E) 

V2 7 = T+(p,q,p = 0,qo = 0)- J fc2(g) g2 (2.90) 

2fj, R 2fj, n 
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In accordance with the perturbation theory methods the corrections of the second order are deter- 
mined by the formula: 



AE H =< V c | AV> 7 | Hf' c >< Hf' c | ^ | y' c > + 
^ < *' c | AF 7 | V nC >< < c | AT> 7 | *' c : 



(2.91) 



with using the Pauli-like eigcnfunctions that are the solutions of the local quasipotential equation with 
the modified Coulomb potential: 



1 I(^c) 2 (l — ^ — ) 

2 mira2 



X1X2, 



*c(0) 



(2.92) 
(2.93) 



(Xi,2 are Pauli two-component spinors), 

b 2 (E)+E 1 E 2 



Veff 



b 2 (E) = mim 2 ^ + (Z«) 2 ^_ 1/2 



-En 



(2.94) 



The corrections of the order (Za) 2 Ep and ^(Za) 2 i?F to the hyperfine structure of a muonium, 
obtained in this approach, are shown in Section IV of the present review. 



3 Calculation techniques for finding the energy 
spectra in the different orders in a 

As one can see from the above-said, a calculation of the fine and hyperfine splitting of the 
energy levels is reduced to finding out the matrix elements of the quasipotential V. In the nonrelativistic 
case, the WF of the system lightly bound can be approximated by the Dirac 8- function. The use of 
the Coulomb IS- state WF gives the opportunity of considering the matrix elements when the non-zero 
momenta p,q^0. However, the relativistic corrections are taken into account more accurately by means 



of the WF in the form of (2.73) when describing bound states. 

In perturbation theory developed in Section II the quasipotential and the interaction kernel include 
the Coulomb Green function essentially. The methods of using this function are different in the problems 
of fine structure and hyperfine structure of the hydrogen-like atoms. In the first case, the main interest 
is the low-frequency region of virtual momenta where an interaction is nonrelativistic. Here, to include 
binding effects in the virtual states of the interaction kernel it is important to consider the block of 
the Coulomb exchanges as a whole, e. g., by means of the expression for the Green function of the 
nonrelativistic Schrodinger equation with the Coulomb potential. 

In analysing the HFS of the energy levels, one should have a possibility to consider contributions 
of the one- or two-transverse-photon exchanges and arbitrary number of the Coulomb exchanges. It is 
sufficient to apply the usual jexpansion of the Coulomb Green function restricting ourselves to the needed 
number of expansion terms Q. 

Up to the accuracy 0(a 5 ), it is sufficient to use the approximate WF 

4%*™® = (2tt) 3 5(p) I cf> c (r = 0) |, E ~ mi + m 2) (3.1) 

when the calculations of the energy of the IS- level, based on the quasipotential V constructed from the 
diagrams of the order a 2 and higher. 



6 One can restrict oneself to the diagrams of the three-photon exchanges in calculating the HFS of the 
ground state in a positronium and in a muonium up to the accuracy a 6 log a. 
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The above-said statement is based on the following fact. The value of squared modulus of the WF, 
which is in the matrix element, has the order a 3 in the coordinate space when r = 0, 

I Mr = 0) | 2 = ~ 0(a 3 ). (3.2) 

Using the well-known representation of the 5 - function: 

1 a 

S(x) = lim a ^ 5— — j, (3.3) 

7T or + x z 

we can find out 

ttS(x) a , . 

-£T = hm ^ {a2 + x2)2 (3.4) 

Thus, 

Snafi | 4> c (0) | 7r5(p) 



lim a ^ 4>c(p) = hm a ^o ,-. 2 2 2 ^ 2 = 87r | C (O) I ■ (3.5) 



Using the formula <5(p) = , which is valid in the case of the spherical symmetry, we get finally 

Um a ^4> c {p) = (2tt) 3 I C (O) | S(p) (3.6) 



It is clear from the form of the Coulomb WF (2.68) that the main contribution in the splitting of the 
energy levels gives the momentum region p 2 ~ Z 2 a 2 fi 2 . As a result, the integrand expansion over p/m 
is equivalent to the integral expansion as a whole over a provided that the integral is finite. 

This property proved to be highly useful in calculations of the terms up to a 5 to the HFS of the 
ground state in a positronium from the diagrams |30[: 



In the calculations of the matrix elements the upper (" big" ) components of bispinors survived and 
it became possible to equate E 2 — m 2 , p, q — in the interaction amplitude corresponding to these 
diagrams. In the calculations of the higher orders, the situation changes because of the singular behavior 
of the integrand expression in the small momentum region. 

Taking the one-transverse-exchange diagram as an example let us consider the extraction of the 
contributions of the order a 2 loga to the Fermi energy of the hyperfine splitting Ep — | ^ ^ < >■ 
The expression of this correction has the form: 



AE^ S =< $' c I F^ 1 



GqKtGq 



F I $c >, (3.7) 



Kt = _ 4TOl M fc ) , (3.8) 
(k 2 - k 2 + it) 



is a kernel corresponding to the one-transverse-photon exchange diagram, 



p ,j\ -. -, (71*0(72*0 ,„ n , 

ri 2 (fc) = 7172 =- — , (3.9) 

K 
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and & c , th e WF in the second order of perturbation theory, can be substituted for <f>c{p) defined by 
Eq. ( |2.68 ) in the calculations up to precision we need, with 4>c{p) being the exact solution of the 
non-relativistic Schrodinger equation with the Coulomb potential for the IS- state. 
The analytical expression of the quantity (3.7) is given by the equation: 

4«V | <mo) | 2 



(27T) 



dpdqdp dqo 
(p 2 + a 2 fi 2 ) 2 (q 2 + a 2 fjpy 



F- 1 (p)F- 1 (q)x 



x / dk dk' [S 1 (pi)S 2 (p2)K T {ko, k' ;p, q)S 1 (q 1 )S 2 (q2)} + S(Po - k )6(k' - q ). 



(3.10) 



Using the Fourier representation of the (5-function and the residues theory, we find: 



AE^ S = 



A%o?\i 2 



(0) 



(2tt) 3 



dpdq ul(p^ul(~p)d 1 d2U 1 {q)u 2 {-q) 
\P~ 9*1 (p 2 + oi 2 iJL 2 ) 2 {q 2 + a 2 [i 2 ) 2 



dte 



i\p-q\\t\ ^(^ e -i\t\(ei g +e 2j> -E-ie) 



■i\t\(e lp +H2q — E-ie) 



(3.11) 



after integration over the variables po, qo, ko, k' . 

Upon extracting the spin-spin interaction from the spin-structure (see the nominator of Eq. (3.11) 
and using a symmetry of spin-structure expressions with respect to the substitutions pi Qj, e. g., 

J dpdq(oiP)(°2q)f{p\ q 2 , (p~ q) 2 ) = ~(ffiff 2 ) J dpdq{pq) f \p 2 , q 2 , (p- q) 2 ), (3.12) 

the integral expression of AE T can be written (the case of unequal masses mj ^ ?ti 2 ): 



AE^ S 



a 3 /z 2 
3tt 3 



I <M0) | 2 < cti(7 2 > / dpdq 



(p 2 + a 2 fi 2 ) 2 (q 2 + a 2 /! 2 ) 2 



(I P — q I + £ ip + ( 2q - E - ie 



Z<p 2 (q 2 ~p 2 ) 



M. 



2q 



+ q 2 (p 2 -q 2 ) 

if- q) 2 



M, 



■2l> 



f-2p + £2<7 £lp + £lg 

(P 2 ~q 2 ) 2 



£2p + £2q £lp + £lg 

(p- q) 2 [M+M+ + M+M+] 



(p-q) 2 {eip + e lq )(e 2p + e 2q ) 
The following notation has been used above: 



(3.13) 



e lp e2 P ei q e2 q M+M+ p M+M+ 



(3.14) 



M tp = £ ip + m ii M iq = e iq + m i- 

As a rule, it is possible to estimate the order in a of each of integrals to the final value of the hyperfine 
shift before an integration. In Ref. |nj it has been shown that the logarithmic correction of the order 
a 2 loga to the Fermi energy appears from the integral 



I st (loga) 



8tt 2 



dp 



dq 



ei P e 2p (p 2 + a 2 fi 2 ) J (q 2 + a 2 fj, 2 ) (p- q) 



p ■ dp 



q-dq p + q 

jlog-. 



£i P £2p{p 2 + a 2 [i 2 ) J q 2 + a 2 [i 2 \p — q\ 2mim 2 
u 

a 6 I st ~ a e log a, 
which is agreed to be called the " standard integral" . 



log a 1 + 0(a), 



(3.15) 
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Additional powers of p or q in the nominator and the additional powers of the factors (p 2 + a 2 /i 2 ) or 
(q 2 + a 2 ii 2 ) in the denominator of the integrand expression lead to the contributions of the order a 6 or 
a 4 , a 5 , respectively. 

The terms in the curly brackets, containing the multiplication p 2 ■ q 2 , lead to the "standard" integral 



in (3.13). Moreover, the difference of the factor (| p — q | +e lp + e 2g — E) from | p — q | in the denominator 
of (3.13) proved to be essential. The main contribution of this term has the order a 4 , but the next order 
terms in the denominator expansion result in the "standard" integral. 

The other diagrams also give the contributions of the order ~ a e log a to the HFS in a muonium. It 
is calculated analogously. The results of calculations, Ref. are shown in Table I. 

Table I. The contributions to the HFS in a muonium from the diagrams of Fig. 2 



No. 


Contribution K l to A£f ; f 


No. 


Contribution K. t to A£f f f 


a 


1/4 


e 


3M 


b 


M + 2 


f 


-2(M + 2) 


c 


9/2 


g 


5/4 


d 


-(M + 2) 


h 


-M 






Total 


2 



In the method based on the amplitude T, when the constituents are on the mass shell pi — 
e ipjP2 = £ 2pi9i = £ i<?j<72 = e 2<ji the problem of correct allowance for the retarding effects appears 
even at the stage of calculation of contributions from the one-transverse-photon exchange diagrams. 
Depending on the way of representation uj 2 in the denominator of a photon propagator 

D « = -^F*-f»- (3 - 16) 

the contribution of the order ~ a 6 log a is different (see Table II). 



Table II. The contributions to the HFS in a two-fermion system from the one-photon exchange 
diagram (lu 2 is zero component of the photon 4 ~ momentum). 





AE^ FS (a (i loga),m 1 ^ m 2 


!S.E^ h b (a 6 log a) , mi = rri2 





E F ^ Mloga- 1 

1 mi mo. z 


\E F o?log a~ x 


[ei P - ei 9 ] 2 


E F ^(M 2 m *)loga- 1 

1 mi mo y mi ' ° 





r i 2 

[ e 2p ~ ^2q_ 


E F ^ a (M 2 m ')loga- 1 

r mi mo, \ mo ' & 





[(eip - ei q )(e 2q - e 2p )] 


E F ^ (M + 2)loga- 1 

1 mi mo. v / j 


E F a 2 log a~ l 



Comparing the result shown in Table II with the result obtained by the first method , the two-time 
Green function method, 

AE* fs (a 6 loga) = E F ^—(^ + — + 2)loga-\ (3.17) 

777,1777,2 TO2 TOl 

we are convinced that it is preferable to use the symmetric form of w 2 in the quasipotential continued 
analytically out of the energy shell (|p| 2 ^| <f| 2 )- 
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—vc 



vc<8> F 



2«) 8K C ® f 



9 h 
Fig. 2. The diagrams contributing to the order of ~ a 6 loga to the HFS in a muonium. 



In the higher order diagrams the situation is more complicated when the above-mentioned method 
of passing to the mass shell is used. Firstly, the way of symmetrization, that is what method should be 
used when going away from the energy shell, is not clear. Secondly, the problem of existence of additional 
singularities complicates essentially the calculations. A loop momentum integration, for instance, should 
be considered in the sense of the main value. The use of various expansions of integrand and the change 
of integration order are problematic in a situation like that. 

Thus, the most correct allowance for retarding effects is the use of the first method for building 
the quasipotential. However, in this case the problem of existence of anomalously large contributions 
~ a 5 log a to the HFS of the ground state in a two-fermion system appears already at the stage of the 
one-photon exchange diagram. In the second method, this trouble does not occur. This problem is not a 
specific feature of the quasipotential approach but it is general for the relativistic bound state theory Q . 

The diagrams contributing to the order ~ a log a to the hyperfine splitting of the ground state in a 
positronium in the direct channel are shown in Fig. 3 . The corresponding quasipotential is presented in 
Refs. [Ill Ejl. 
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These contributions have already been mentioned to be caused by the infrared behavior of matrix 
elements of the quasipotential. The existence of the iteration terms for each of the reducible diagrams 
improves its behavior in the infrared region. This fact permits one to avoid summing the ladder diagrams 
in any selected order of a. The parameter a/j, is found out to behave as the regularization factor of the 
infrared singularities when the calculations with the exact Coulomb WF (2.68) are carried out in the first 
variant of the quasipotential approach. The cancellation of these anomalous terms [[l5] is displayed in 
the Table III. 



T 



a b c 




e f 



-2® v c ® F 

L 

g 



F 



h 

Fig. 3. The diagrams considered in analysing the anomalous contributions of the order ~ ct'loga. 
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Table III. The cancellation of the anomalous contributions of the order ot'log a to the HFS 
in a positronium. 



No. 


AE(a 5 loga) 


No. 


AE(a 5 loga) 


a 





e 





b 


'fE F loga 


f 


-fE F loga 


c 


-7^E F loga 


g 





d 


^E F log a 


h 









Total 






The new corrections to the hyperfine splitting of energy levels in a muonium have been obtained by 
Eides et al. |3J| — j35|. In these articles, in particularly, the corrections of the orders ~ a(Za)^- and 
~ Z 2 a(Za)^- to the Fermi energy have been calculated by means of the method of the effective Dirac 
equation (EDE)i. Let us consider the selection of these contributions from the diagrams with radiative 
photons. The remarkable feature of these articles is the use of the Fried- Yennie gauge |3j|-[^lJ for the 
photon propagator 

2V = -f^(^ + 2-gH (3.18) 
q z + le + it 

The infrared singularities are softened in this gauge. Any diagram with radiative corrections has a 
softer behavior near mass shell than the corresponding "skeleton" diagram. An attractive property of the 
Fried- Yennie gauge is the possibility of carrying out the renormalization procedure on mass shell without 
introduction of the unphysical photon mass A. This feature makes it easy to estimate integrals appearing 
in the problems of energy levels of the hydrogen-like atoms. 

In the Fried- Yennie gauge the renormalization constant for the WF, Z2, is infrared finite and the 
renormalized self-energy operator has a soft behavior on mass shell 

£$(p) = (P-™) 2 (~|^)(1 + 0(P)). (3-19) 

2 2 
m — p 

P = 



4%) = ip-m)- 

_A 

rn 



This is a distinctive feature of the mentioned gauge from, e. g., the Feynman gauge 

<— -logp+1 , (3.20) 
vi l m 

< p< 1. 

As for the vertex function, the term corresponding to the fermion anomalous magnetic moment: 

- f <V#- (3-21) 
2tt 2m 

has a most hard behavior. However, redefining the renormalized vertex operator by means of 

A>i, P2 ) = 7 M A(0,0) - + AW( Pl ,p 2 ) (3.22) 



7 Recently, these authors have calculated the corrections of the order a 2 (Za)EF ■ See Section IV and 
l-M for the details. 
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one gets the expression: 




3a (p — m)p 



(3.23) 



47r m 2 



which is valid when the transferred momentum is zero and /) <C 1. This expression agrees with the 
self-energy operator asymptotics owing to the Word identity. The contributions of the term according to 
the anomalous magnetic moment is analyzed separately. 

Let us trace the selection of graphs for the calculation of the corrections of the order ~ a(Za)^-Ep- 
Five diagrams: 



exhaust the contributions to the EDE kernel, connected with the mass operator. 
Here, 



with 5*0 being the free particle propagator; g = A 2 • S being the projector onto the muon mass shell 
multiplied by an electron propagator. 




s 



A 2 -S 
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The simplest diagrams with the vertex correction: 



also give contributions of the order ~ q(Z«)^. 

Moreover, there are diagrams with spanned many emitted photons: 



and the diagrams of the second order of perturbation theory: 



21 



In analysing the graphs, entering into the EDE kernel, with the aim of finding the corrections of the 
order ~ a(Za)^- to the Fermi energy, it has been determined that the contributions of this order come 

from the diagrams of the gauge-invariant set only (see Fig. ^)Q. Moreover, it turned out to be possible 
to restrict oneself by the approximated WF ( |3.l| ) in the matrix element. In other words, the matrix 
elements are to be calculated with taking into account the upper ("big") components of electron and 
muon spinors, neglecting momenta of wave functions inside the diagrams. These conditions were named 
the "standard conditions" by the authors of p3[-p5[. 



Fig. 4. The complete gauge-invariant set of the diagrams for the calculation of the recoil corrections 
of the order ~ a{Za)Ep and ~ Z 2 a(Za)Ep to the Fermi energy of the HFS in a muonium. 

8 The diagrams, where the radiative photon spans more than two exchange photons, don't contribute 
to the terms of the order ~ a(Za)Ep in the FY gauge. 
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The complete gauge-invariant set, presented by Fig. 4i leads to the infrared and ultraviolet finite 
matrix element in a sum. Therefore, further calculations can be done using any convenient gauge both 
for the exchange photons and for the radiative photons. It has been deduced by direct calculations that 
the anomalous magnetic moment does not lead to the corrections of the order needed. It is considered 
to be subtracted from the vertex operator. The iteration diagrams of the EDE method, marked by x 
on the muon line, are analogous to the diagrams of the quasipotential approach in some sense. Like 
in the quasipotential approach there are two iteration diagrams for the two-photon exchange diagrams 
with the radiative insertion (t^Ft^ and t'^Ft' 4 '), there are "non-compensated" iteration diagrams 
in the EDE method. They are generated by the BS kernel including the inverse muon propagator. The 
iteration diagrams turned out , Ref. to be canceled by the vertex correction diagram in the order 
under consideration after building the EDE kernel and perturbation theory for finding the energy levels. 

As a result, the following expression is obtained for the contribution from the diagrams of Fig. 4a: 



87r 2 /t 2 
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= SE^ +SEz2, (3.24) 



x and y are the Feynman parameters, fj, = 

1+3 



hi(x,y) 



y 



al(x,y) 



x 2 + X 2 
(1 - x)y' 



2(1 



x 2 + X 2 



(3.25) 



A is a non-dimensional infrared mass of a radiative photon in the electron mass units. 

The main contribution of the order 1/ p, to the integral comes from the residue in the muon pole, what 
corresponds to the muon motion on the mass shell. The leading infrared singularity, which is proportional 
to A -1 / 2 , is also connected with this residue; the other terms are logarithmic divergent only. Therefore, 
it is convenient to separate the calculation of the on-shcll contributions and the rest of them. 



5E s (m. s.) = — (-2I X + 
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where 



h = ^ldx 



/ \ 3 / 2 

(^) a-)" 2 



A-V 2 



(3.26) 



(3.27) 



is the infrared divergent integral which is subtracted after summing the singular contributions of the 
diagrams Fig. J^a, 4b, 4 C - 

By means of a number of mathematical contrivances (see [}5| ) like dividing the integration region 
into two parts, small and large momenta, after the integration over the angular variables, the subtraction 
of the pole contribution in the integrand and use of various identities, the expressions for the contributions 
from all the diagrams of the gauge-invariant set have been obtained, Fig. 4'- 
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Another method, which gives the opportunity to calculate the logarithmic contributions of the orders 
a 3 — log 3 — to the HFS in a muonium, is the renormalization group method, Ref. 

The contribution of the radiative-recoil corrections is equal to AE = —EpRfn where the quantity 
is calculated in the lowest approximation from the diagrams of two-photon exchange. In Ref. 0], the 
diagrams with the radiative insertions into an electron line, photon line and to an electron vertex have 
been considered. The following logarithmic contributions are known from the above paper: 



R)P = -log— + {-? — 

h IT TO,, TO,, 7T TO,, 



2 TO e 31 TO e 

2log h — log 

to m 12 



(3.29) 



which arise owing to momentum integration in the region m e « r < mt, that is in the asymptotic 
region for the contribution of the electron vacuum polarization to the photon propagator. The muon 
loop does not give contribution under momentum integration in the region k 2 < m 2 ^. The estimation of 
the term of the higher order of perturbation theory is given in Ref. |Q . 

If we consider some physical quantity R, calculated by perturbation theory, the following condition 
is to be fulfilled: 

£ = 0, (3.30) 

where the variable r = —(3$logj^ characterizes the used renormalization scheme (RS), fi is an arbitrary 
parameter of mass dimension, A is the scale parameter [f45{ , /3 = 2/3 is the first coefficient in the 
renormalization group equation with the running coupling constant g, 

A^ = /?(5) = /3o.9 2 +/3i5 3 + - (3.31) 
Accordingly, for the quantity R in the second order of perturbation theory which is written in the form: 

=r g(l + r l9 ), (3.32) 

we have: 

dR^ 

^=0(g 3 ). (3.33) 

It follows from here that 

Thus, we can see that tq does not depend on the selection of the RS and T\ = t + pi, where the 
constant p\ can be calculated provided that r% is known for some of the RS. 

The dependence of the next (third) term of perturbation theory on r is to be arranged so as to 
compensate the dependence R^ 2 > on r up to the order g 4 , 

d(R& + f)( V) _ 4 



dr 

Consequently, 



0(5 4 )- (3.35) 



^-=r (2r 1 + ^). (3.36) 
or p 



After integration using Eqs. (3.34) we get 



^ (2) =r ri(ri + §-) + co7isi. (3.37) 

Po 



9 The mentioned method has also been used in the calculations of terms depending on the logarithm 
of the mass ratio to the value of the anomalous magnetic moment |fl3| . 
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Then, to define an arbitrary integration constant it is necessary to set the "optimal" RS in which the 
quantity is the closest to R. In Ref. p2| it was given by the condition 



n^( ri ,r) \ T=Top = 0,r opt = -Polog(^) 



If we use the RS on the mass shell (ri = Ki) as the initial RS, then we have 

r° pt = n(r pt)- 



Po Po 



(3.38) 



(3.39) 



We obtain from ( 3. 29] ) for a physical quantity AE, the energy of the HFS in a muonium, the following 
expression: 

■ > .,. • > i 

(3.40) 
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Also, from (3.39) we have 
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As a result 
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(3.41) 
(3.42) 

(3.43) 



The first term in (3.43) is equal to the result obtained in |33|] by a direct calculation. The value of A in 
(3.41) depends on the selection of the "optimal" scheme. Its contribution to the HFS can be estimated 
only as < — AE < 1 kHz. 

The remarkable method of calculation of the logarithmic in a corrections to the spectra of the QED 
systems has been proposed in [^J. The logarithmic in a corrections are pointed out in these papers 
to appear from the logarithmic divergent integral, with the contribution being given by the momentum 
region 

fia<q<n (3.44) 

(fx is the reduced mass, fi = ™+ l M ) ■ The lower limit is the characteristic momentum for QED bound 
states (see the discussion of the WF ( 2.68 ), the upper value corresponds to the limit of applicability of 
the nonrelativistic approximation. 

A shift of the level with the quantum numbers n, I has been calculated in |47j| as the matrix element 
of the operator V(q): 



V^(q) = [A + B(a 1 a 2 )] 



N 

na /J, 
log- 



mM 



9 



(3.45) 



with the WF's which are similar to the ones presented by Eq. (3.1). A calculation of the operator V(q) 
is reduced in this approach to a calculation of the on-shell scattering amplitude from the diagrams of the 
order a , if we are interested in the logarithmic corrections only. 

The authors of the paper |47{ have used the usual quantum-mechanical perturbation theory with 
the substitutions 

| i >< i | 



i 

in the intermediate states of an interaction operator. 



(3.46) 
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The coefficient of the term log — has been overestimated in 
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is: 




J (2tt 3 )J (2tt 3 ) k 2 k' 2 (k-k'y (E - eik - e 2k )(E - e lk . -em) ' 



With the results from other diagrams of the Coulomb photons exchange this leads to the level shift 



The theoretical values of the decay width of o — Ps and p — Ps and of the energy levels of fine and 
hyperfine structure have been obtained by I. B. Khriplovich et al. Their results are discussed in the 
next section. They are compared with Fell's result, Ref. jl8| calculated on the basis of the relativistic 
two-particle equations. 
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Abstract 

The present review includes the description of theoretical methods for the inves- 
tigations of the spectra of hydrogen-like systems. Various versions of the quasipo- 
tential approach and the method of the effective Dirac equation are considered. 
The new methods, which have been developed in the eighties, are described. These 
are the method for the investigation of the spectra by means of the quasipotential 
equation with the relativistic reduced mass and the method for a selection of the 
logarithmic corrections by means of the renormalization group equation. The spe- 
cial attention is given to the construction of a perturbation theory and the selection 
of graphs, whereof the contributions of different orders of a, the fine structure con- 
stant, to the energy of the fine and hyperfine splitting in a positronium, a muonium 
and a hydrogen atom could be calculated. 

In the second part of this article the comparison of the experimental results 
and the theoretical results concerning the wide range of topics is produced. They 
are the fine and hyperfine splitting in the hydrogenic systems, the Lamb shift and 
the anomalous magnetic moments of an electron and a muon. Also, the problem 
of the precision determination of a numerical value of the fine structure constant, 
connected with the above topics, is discussed. 



* Submitted to "Physics of Particles and Nuclei" on the invitation, 
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4 Comparison of theoretical and experimental 
results 



4.1 Decay rate of a positronium 

Quantum-electrodynamic systems, consisting of a particle and an anti-particle, have specific features. 
Apart from a scattering channel an annihilation channel appears in this case. A positronium atom, which 
is a specimen of these systems, has no stability. The life time of a positronium (or the decay rate) is 
the subject of precise experimental and theoretical investigations. The charge parity of a positronium, 
C = ( — l) i+s (L is the eigenvalue of an angular momentum operator, S is the eigenvalue of a total 
spin operator for the system under consideration), is a motion constant. Consequently, all its states are 
separated into the charge-even states (S — 1) and the charge -odd states (S = —1). The positronium 
total spin is also conserved and the energy levels are classified as singlet levels (S = 0, a parapositronium) 
and triplet levels (S = 1, an orthopositronium). The S- state (L = 0) parapositronium has a positive 
parity and the S- state orthopositronium has a negative parity. As a consequence of conservation of a 
charge parity in electromagnetic interactions a parapositronium disintegrates into the even number of 
photons and an orthopositronium into the odd ones. 

At present, essential disagreement exists between the theoretical and experimental values for the 
decay rate of an orthopositronium. The theoretical predictions are Q-[E2| 



Tf eor {o-Ps) 



= To 



a Q mc 2 2(ir 2 - 9) 
h 9?r 
-1.984(2) \ 
-1.9869(6)J 



7.038 31(5) [is- 1 , 



I -An- - 2l - 1 1 



-log a 



-uHoga- 1 + B 3 (-) 2 + 



4 2 
-9 C(2)+ 3 + — V 



(4.1) 



where 



.4 
A\ 



-10.266 ± 0.011, 
-10.282 ±0.003. 



The last experimental mcasurings are |33[ [54)0 

r||| (o-Ps) = 7.0514(14) [is^ 1 
rgjt (o-Ps) = 7.0482(16) us' 1 . 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



The result of Ref. [g3| has 9.4 standard deviation from the predicted theoretical decay rate and the result 
of Ref. |Q has 6.2 standard deviation. The coefficient = 1 in 0(a 8 ) term can contribute 3.5 • 10 _5 //s _1 
(or 5ppm of T^) only. To remove the above disagreement the coefficient B% should be equal to about 
~ 250 ± 40. It is very unlikely, indeed, but this opportunity has been pointed, out in fl54| and cannot be 
rejected a priori. The calculation of the B% coefficient is very desirable nowlij. 

For the first time, the main contribution in the orthopositronium decay rate has been calculated in 
Ref. pf: 

2 



T (o-Ps) = -2Im(AE 37 ) = 

y7T 



9)ma b = 7.211 17 /xs" 



(4.6) 



The corrections of the O(a) order to this quantity have been calculated in a numerical way ||50| , p2| , |56| , 
B7[ at first, but later some of them have been found analytically, Refs. to, M,[M-|6fj], by using the 



11 See Table IV for the previous experimental results. 

12 Some estimations of the corrections of this order have been done in Refs. 
AE= 28.8(2)(a/7r) 2 r . 
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Feynman gauge. The corrections arising from the diagrams with self-energy and vertex insertions have 
been calculated by Adkins 59, [30| 
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where 



and 



/? = / dx 1 ^- — — [C(2) - Li 2 {l - 2x)} = -1.743 033 833 7(3), 
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1 °° 1 



n 4 2 r ' 



0.517479 061674, 



C(2) = — , C(3) = 1.2020569032, 
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D = 



1 



2-w 



- 7b + log (An) 



(4.10) 

(4.11) 
(4.12) 
(4.13) 



is the standard expression of a dimensional regularization (2w is the space dimension.). The above results 
is co-ordinated with Stroscio's result IB8I when 



r - [-D - 4 - 2log(\ 2 /m 2 



(4.14) 



is added to the last one. This is necessary to do because of different regularization procedures used in 
§§ and gg, respectively. 

Recently, calculations of these corrections have been completed, Ref. pU, in the Fried - Yennie gauge 
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(4.15) 



(4.16) 



The contributions from the remained diagrams, with a radiative insertion in a vertex of an internal 
photon; with two vertices spanned by radiative photon; the diagram taking into account binding effects 



2 



and the annihilation diagram (see Fig. 1 in |n|-b), have been calculated numerically. As a sum the 
0(a) corrections are jointed to give 

7 

in rv rv 

— i- [-1.987 84(11)] = r Q - [-10.2866(6)] . (4.17) 

ThenEl 

r^^(o-Ps) = 7.038 236(10) ^s" 1 . (4.18) 

The above result is the most precise theoretical result at present. 

To solve the existing disagreement between theory and experiment, the 5- photon_mode of o — Ps 
decay and the 4- photon mode of p— Ps decay have been under consideration in (6^, ^|tl The following 
theoretical evaluations were obtained: 



(o - Ps) 



T 3 {0-Ps) ' TT 
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0.177(-) 2 - 0.96- 10~ 6 , (4.19) 



and 



F f, } P = 0.274(-) 2 ~ 1.48- 10~ 6 , (4.20) 

L2(P—PS) TT 



rP(o-Fs) = 0.018 9(ll)a 2 r , (4.21) 
rfi (p-Ps) = 0.013 89(6)ma 7 . (4.22) 

They are in agreement with one another and with the results of the previous papers [|64|Ei 

r 4 @(p-Ps) = 0.013 52 ma 7 = 11.57 ■ 10~ 3 s~ 1 . (4.23) 

In the connection of the present situation concerning the decay rate, investigations of alternative 
decay modes for this system (e. g., o — Ps — > 7 + a, a is an axion, a pseudo-scalar particle with mass 
m a < 2m e ) are of present interest, Refs. ^FJ-[£l). In the paper |6^], the following experimental limits of 
the branching of the decay width have been obtained: 

Br = T{ ° P f,^ 7 + Q) < 5 ■ 10- 6 - 1 • 10- 6 (SOppm), (4.24) 
I (o — Ps — > 37) 

provided that m a is in the range 100 - 900 keV. In the,-£ase of the axion mass less than 100 keV (This is 
implied by Samuel's hypothesis |7(J. According to it E3 m a < 5.7 keV, g ae + e - ~ 2 • 10 -8 ) the limits of 
Br are the following jnj : 

Br = 7.6 -10" 6 , if m a ~l00keV, (4.25) 
Br = 6.4 ■ 10~ 5 , if m a <30keV. (4.26) 

These limits are about 2 orders of magnitude less than the value which is necessary to remove the 
disagreement. 



13 The uncalculated yet 0(a s ) corrections are not taken into account here. 

14 As a consequence of conservation of an angular momentum and isotropic properties of the coordinate 
space, an orthopositronium has to decay into the odd number of photons and a parapositronium into the 
even ones, see above. 

15 The result of Ref. |65[] is not correct, four times less than the above cited results. The explanation of 



this was given in 63 . 

16 The proposed values don't cause the contradictions in comparing theoretical and experimental results 
of the anomalous magnetic moment (AMM) of an electron. 
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Finally, a decay o — Ps — > nothing (that is into weak-interacting non-detected particles )0 has been 
investigated in Ref. [173). The obtained result 



r(o — Ps — > nothing) 
r(o- Ps 3 7 ) 



<5.8-10~ 4 (350 ppm) 



(4.27) 



excludes an opportunity that this decay mode is the cause of disagreement between theory and experiment. 

The decay of o — Ps into two photons which breaks the CP- invariance, as has been mentioned in 
|74j, f75|, was experimentally rejected in |76| E3. 

It should be marked that the contribution of a weak interaction has been studied in j78j. However, 
because of the factor rr^jM'k, it cannot influence the final results. In the cited articles the weak decay 
modes have been estimated as 

r(p - Ps 



37) r( - Ps -> 47) 



a(G F m 2 e g v ) 2 



r(p - Ps -> 27) r(o - Ps -> 37) 

where is the Fermi constant for a weak interaction, 

SV = 1 - 4sin 2 e w ^ 0.08, 
6(y is the Weinberg angle. The present experimental limits are [ f79| , 

r(p - Ps -► 37) 
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-27 



r(p- 
r( - 



Ps 
Ps 



27) 
47) 



< 2.8-10" 



< 



• 10" 



r(o- Ps -> 3 7 ) 

In Table IV all experimental results for the o — Ps decay rate , known to us, are presented 0. 
Table IV. 



(4.29) 

(4.30) 

(4.31) 
(4.32) 



Year 



1968 
1973 
1973 
1976 
1976 
1978 
1978 
1978 
1978 
1982 
1987 
1987 
1989 
1990 



Reference 



[84] 
[85] 



[8£] 
[9C] 
[91] 
[92] 
[93] 
[94] 
[95] 
[53] 
[54] 



Ps) 



flS 



7.262(15) 

7.262(15) 

7.275(15) 

7.104(6) 

7.09(2) 

7.056(7) 

7.045(6) 

7.050(13) 

7.122(12) 

7.051(5) 

7.031(7) 

7.0516(13) 

7.0514(14) 

7.0482(16) 



Error, ppm 



2070 
2070 
2060 
840 
2820 
990 
850 
1840 
1680 
710 
1000 
180 
200 
230 



Technique 



gas 
gas 
gas 

powder SiC>2 
vacuum 

gas 

gas 
vacuum 
vacuum 

gas 
vacuum 

gas 

gas 
vacuum 



Regarding the results for the decay rate of a parapositronium, the situation was highly favorable until 

17 Like that Glashow (72) spoke out the hypothesis of the decay into invisible "mirror" particles. 
18 The physics ground of these speculations is possible existence of an unisotropic vector field with 
non-zero vacuum expectation pq , with which an electron and a positron could interact 

C = gTpO al3 ipA a nP, (4.28) 

C is the interaction Lagrangian. 

19 The results of the papers |M[ ^2| and p3] could be accounted as rough estimations. 
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the last time. The theoretical value, found out as early as in the fifties, Refs. pa, B7J, is 



r theor {p _ ps) = _2Im(AE 2l ) = - 



1 a 5 mc 2 



a , 7r . 

7T 4 



= 7.9852 ns" 



(4.33) 



2 ?i 

The above value, confirmed in |9^, |9^], coincides with the direct experimental result with good accuracy 

Tjg (p - Ps) = 7.994 ±0.011 tis" 1 . (4.34) 

The experimental values of the parapositronium decay rate are shown in Table V El 

Table V. 



Year 


Reference 


T 2 (p — Ps), ns 1 


Error, % 


Technique 


1952 
1954 
1970 
1982 


I 
I 
I 


100 
101 
102 
93j 


1 
1 
1 


7.63(1.02) 
9.45(1.41) 
7.99(11) 
7.994(11) 


13 
15 
1.38 
0.14 


gas 
gas 
gas 
gas 



In Refs. |5lL it has been pointed out that it is necessary to add the logarithmic corrections in a 
to Harris and Brown result. In the paper [[|7|-a these corrections to 1^3(0 — Ps) and T 2 (p — Ps) have 
been re-calculated, with the result of the decay rate of a parapositronium differing from the one found 
out before, Refs, pi |9l: 



rpi(p- Ps,a 2 loga) 
T^^(p-Ps,a 2 loga) 



ma 
ma 5 



■ 2a 2 log a , 



• —a 2 loqa . 
2 3 y 



(4.35) 
(4.36) 



Finally, we would like to mention a quite unexpected result, presented in Remiddi's (and collabora- 
tors) talk [103]. The calculations carried out by them lead to the additional contribution 



rpl(p- Ps,aloga) = ^(-)2loga, 

2 7T 



(4.37) 



which is explained by them to appear as a result of taking into account the dependence of the interaction 
kernel on the relative momenta (see, e. g., Fig. 1). 

The above-mentioned leads to necessity to continue calculations of the decay rates of an orthopositro- 
nium as jwell as a parapositronium using more accurate relativistic methods, e. g., the quasipotential 
approachEil 



4.2 Hyperfine splitting 

4.2.1 Positronium 

Comparison of theoretical and experimental results of the hyperfine splitting of the ground state 
of a positronium and a muonium was considered for a long time as correction of our understanding 



20 The branching of the decay rates of a para- and an orthopositronium, p^7~p~j ; had been measuring 
in the experiments of 1952 and 1954. The presented results are recalculated by means of the first direct 
experimental value, Ref. @, T 3 (o - Ps) = 7.262(15) ps" 1 . 

21 It also deserves an atte ntion The new approach to the positronium lifetime calculation, proposed by 
A. A. Pivovarov, Ref. [ 104 1 , also deserves some attention. 
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the bound state problem. At first, the quantity of this splitting was estimated in a positronium as 
9.4 • fCT 4 ± f .4 • iCT 4 eV in 1951, Ref. p5|. The last measured result is, Ref. pig: 



AE c h f s (Ps) = 203 389.10 ± 0.74 MHz (3.6 ppm). 



(4.38) 



In Table VI all published values of precise experimental measurings of the hyperfine splitting of the 
ground state of a positronium are presented. 

Table VI. 



Year 



1952 
1954 
1955 
1957 
1970 
1972 
1975 
1975 
1977 
1983 
1984 



Reference 



[107| 



[108| 



[109 



[110| 



[102| 



[1H| 



[112| 



[113| 



[114| 



[115| 



[106 1 



AE, GHz 



203.2(3) 

203.350(50) 

203.380(40) 

203.330(40) 

203.403(12) 

203.396(5) 

203.3870(16) 

203.3849(12) 

203.384(4) 

203.3875(16) 

203.38910(74) 



Error, ppm 



1500 
250 
200 
200 
60 
24 
8 
6 
20 
8 
3.6 



All the conducted experiments are based on the techniques using an observation of Zeeman transitions 
in Ps and the further substitution of the results into the well-known Breit-Rabi equation. AE, the energy 
of the HFS, is deduced from it. 

At present, the theoretical result determined [ 



lLM, [ne]-[i20| is 



A T^theor 



(Ps) 



7 a. 8 log2, 5 2 , _-, 2 \ 



203 400.3 MHz. 



(4.39) 



The coefficient 1 of the term of an order ~ a 6 can contribute « 18, 7 MHz to the energy of the HFS. The 
estimated uncertainty is almost 50 ppm, an order of magnitude greater than the experimental one. After 
calculations of the corrections of an order 0(a 2 ) and 0(a 3 log a) to the Fermi energy the theor etica l erro r 
would decrease to lppm. The work in this direction has been started since the seventies [Q, 121 -|13C|. 



The first contributions of this order have been calculated from-.the diagrams of a one-photon annihi- 
lation with a polarization insertion of the fourth order, Ref. [122p3, 



A£&£ = ia 2 i?oo(-) 2 

Z 7T 



324 8 SU 4 y 



35^ 2 



32 



= -2.78 MHz. 



(4.40) 



The contribution of a three-p hoto n annihilation [123], which has also been calculated analytically 
has recently been corrected, Ref. [ 124 1 , 



AE. 



hfs 
37 



a R a 



r C(3)--C(2)io 5 2- 



|C(2) - 4/o 5 2 + | 



;C(2) 



(4.41) 



-0.969 MHz. It has been confirmed 



This expression contributes numerically the small value Re(AE^) 
by the authors of p| in Ref. p5 |. 

The analytical expression of the contribution from the two-photon annihilation diagram is 



AE^ S = 



2tt 2 



35 



,41 



85 



1 + ^ + (t" + * ) l °9 2 ~ TC(3) ~ i<5 ~ -r) 



(4.42) 



22 The only contribution coming from electron-positron loops is essential in the case of positronium 
because of the smallness of m P . 
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At present, it gives the greatest contribution Re(AE2 1 ) = —13.13 MHz numerically as compared to other 
corrections of this order. 

It is to be mentioned that in calculating corrections of the order ~ a 6 different authors used different 
approximate methods and, consequently, the comparison of their results is rather difficult. So, the authors 
of Ref. [128 1 investigated the case of the static interaction kernel, which is the fourth component of the 
vector potential independent of the relative times (the method of an effective potential) . They found the 
result up to the third order of perturbation theory 



AEjtt = -4a 4 i? 



where 



12 



-hl + log\a) + -G+-{l-AF) 



a A R c 



24 



log a 



0.031 



G = 



tg 1 p 



dp = 0.915 96. 



(1.9 + 0.3) MHz, 



(4.43) 



(4.44) 



is the Catalan constant; 



F = 



{tg^pf 



dp = 0.3897.. 



(4.45) 



In |44|, |l3Cj the contributions of various diagrams have been calculated by numerical methods. The 
importance of finding the contributions ~ a 6 from the diagrams uncalculated until now was pointed out 
(see Fig. 7 in Ref. @). 



f. [131 1 of the contributions from a 



For the completeness let us mention the recent calculation, 
weak interaction to the HFS of the ground level of positroniumc3 . The obtained result, AE^J e s ak 
4.76 x 10 -14 eV = 1.15 x 10 -5 MHz, is far from the present experimental precision. 



4.2.2 Muonium 



In the previous reviews [J35J, [44|, |135[ -[137 the following theoretical result for the HFS of the ground 
state in a muonium was given: 



AE^(Mu)=E F (l + a^ 



2 1 + a,, 



(4.46) 



with 



^,2 J) ^§^Z,u 



log(Za) — log 4 - 



281 
480 



2:) 



The first calculations have been done some years ago [132 -| 134 
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a(Za) 2 
3{Za) 

IT 



(15.38(29))- 



a 2 (Za) 



D, 



2 lo 9 + (Za) 

% m e (to 



to m ) 2 



-2log{Za)- 



8log2 + 3^ 
18 



(4.47) 



(4.48) 



The classical works 1 138 - |14C| were devoted to calculations of th e non -recoil corrections without tak- 
ing into account of finitencss of mass of a heavy particle. In Refs. [^4, 130 the result for the corrections of 
the order a(Za)Ep has been confirmed and the contribution of the order a(Za) 2 Ep has been obtained nu- 
merically. The leading recoil corrections have been calculated by many authors [^5. 11L 12C, 135 , |141|- 
[143]. In Refs. |m| [35f the contributions of the diagrams with radiative inseijtipns in electron and muon 
lines, which depend on — , the mass ratio, have been calculated analytically 2 -* 



i5 M ( — -; electron line) 



a(Za) m e 



rn 



rn,, 



15 Hlfj 

—log — 

4 to r 



6C(3) + 3Tr 2 log2 



TV 

T 



17 



; muon line) 



Z 2 a(Za) m e 



-C(3)-3^ 2 /o 5 2 



39 



(4.49) 



(4.50) 



In Ref. [144 1, the above-mentioned result has been confirmed for the contributions of an electron line. 
The calculations were carried out in the Fried - Yennie gauge for radiative photons. 



The contributions of the vacuum polarization diagrams have been calculated earlier, Refs. [145, [14 

T 2 



rn 



7r m b 



y y m' 3 yK m' 9 3 



The leading logarithmic corrections with respect to m e /m^ have been considered in |4^, ^[El: 



_ ,m e , m u . a 2 (Za) m e 
8^-llog-±) = V^' — 



\WC^)-\log 2 i^) 

3 TO e 3 TO e 



(4.51) 



(4.52) 



Moreover, the program of calculation of the puxe radiative corrections of the order a 2 (Za)Ep (that 



is of the D coefficient in (4.47) will be finished soonEl Refs. 



The a 2 (Za)Ep corrections, induced 



by the diagrams with the insertions of vacuum polarization loops into external photons (see Figs, la 
lc in the cited paper) have been calculated in |36 -a 

a 2 (Za) 



AE(a 2 {Zct)) 



-Ep 



38 

3 

15 



7T 

916391 
37800 



f 4, 2 l + \/5 20 /- 1 + ^/5 608, n it 2 

-VZi og — + —i og2 + T 



n a 2 {Za) „ 
-2.23 — -Ep 



-1.2 kHz JorMu 
-0.34kHz forH 



(4.53) 



induced by the diagrams with the insertions of vacuum polarization loops into radiative photons (see Fig. 
1 in the cited paper) have been calculated in [^6|-b, 



AE(a 2 (Za)) 



a 2 (Za) 



E, 



149 
270 



2 

9vr 



dq ■ D(q) 



-arctg 



2q 



n 



24 0riginally, these corrections have been found by numerical methods, Ref. 13C ] . 

25 The result of the paper J42|, obtained by the technique of the preceding section , has been completed 



26 Recently, Ref. |3q], Kinoshta has presented the preliminary result of calculations of the last remained 
diagram of the order a 2 (Za) (with cross virtual photons). It is AE(a 2 (Zaj) ~ (— 0.64±0.06) - Ep = 
—0.353(33) kHz for Mu, what gives the opportunity to reduce the theoretical error in ( p~59| ) to 0.17 kHz. 



s 



1 



2927 10169 



1 - q \ 4 1 + q 2400 
0.17 kHz forMv> 
0.054 kHz forH / 



3600 



-0.310 742. 



a 2 {Za) 



(4.54) 



D(q) is the total elliptic integral; 

induced by the diagrams with the insertions of the light-to-light scattering sub-diagrams (see Fig. 1 in 
the cited paper) have been calculated in Ref. |36[]-c, 



AE(a 2 {Za)) ~ -0.48213. 



a 2 {Za) 



E F 



-0.26 kHz forMu\ 
-0.084 kHz forHj 



(4.55) 



Let us mention the first calculated recoil corrections of the second order in m e /m^, Ref. [ |t47[, 
corrections of the (Za) 2 Ep and — (Za) 2 Ep orders, calculated in the quasipotential approach, Ref. J24f 2 



AE = E F < 1 + (Za) 



m M ) 2 



19 

y 



1 2 , 

7T" 

72 3 



the contributions of the hadron vacuum polarization (44|, |148] ] , 



6 U IM§I (hadrons) = 3.7520 ± 0.2373 ~ 0.250 ± 0.016 kHz 

p \7r/ mi 



and the estimations of the weak interaction contributions, Refs. [132, 133 1 

3 



AE(weak int. 



-(Za) 1 GF"m e m tl EF — 0.065 kHz. 



(4.56) 



(4.57) 



(4.58) 



4a/2tt 

As a result of inclusion of the above contributions, the theoretical predictions for the HFS of the 
ground state in a muonium is 

A£^7(Mu) = 4463 303.0(0.2)(1.3)(0.6)/c#z. (4.59) 



The first uncertainty arises from the experimental error in determination of a, Ref. |149], the second one 
is from the experimental error in determination of the ratio of electron mass and muon mass, Ref 15C], 
the third one is from the coefficient D that is not finally calculated. 

Regarding the problem of experimental measuring of AE^ P S (Mu) the published values for the HFS_of 



the ground state in a muonium are summed in Table VII beginning from its discovery in 1960, Ref. [151 Ej 
Table VII 



Year 



1962 
1964 
1964 
1969 
1969 
1970 
1971 
1971 
1972 
1973 
1975 
1977 
1980 
1982 



Reference 



[153| 



[154 1 



[156 1 



[157| 



[1581 



[159| 



[160 1 



[161 



[162 1 



[163| 



IL§4 1 



[165| 



[150| 



AE, kHz 



4461300(2200) 
4463 330(190) 
4463150(60) 
4463 260(40) 
4463 317(21) 
4463 302.2(8.9) 
4463311(12) 
4463 301.17(2.3) 
4463 240(120) 
4463 304.0(1.8) 
4463 302.2(1.4) 
4463 302.35(52) 
4463 302.90(27) 
4463 302.88(16) 



Error, ppm 



493 
43 
13 
9.0 
4.7 
2.0 
2.7 
0.5 
26.9 
0.4 
0.3 
0.12 
0.06 
0.036 



27 The three-photon exchange diagrams have not taken into account there. 

28 The result of Ref. pi of 1961 has the meaning of rough estimation, AE = 2250 - 9000 MHz. 
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We conclude that the agreement between theoretical and experimental results is excellent@. On 
the basis of this fact one can get the value of the fine structure constant (see (4.199) by comparing the 
expression (4.46) and the recent experimental result of Table VII. 



4.2.3 Hydrogen 

The analytical result of calculations of the quantum-electrodynamics non-recoil corrections to the 
HFS in a hydrogen is obviously the same as in a muonium (4.46, 4.47 The numerical value for it is 



AE(QED) = 1420.451 95(14) MHz. 



(4.60) 



The recoil correcti ons and dynami cs correction caused by the nuclear structure h ave been c a lcula ted 
in the classical works [141 ,[167|-[169| which have been completed by the results of |l66| ,[170|-[ [l75| l 3 -il: 



AE '(structure) 
5 p (Zemach) 
Sp(recoil) 
dp(polarizability) \ 



= Ef [5 P (Z emach) + S p (recoil) + S p (polarizability)] . 

= -2^a < r p >~ -38.72(56) ppm, 

~ 5.6&ppm, 

< 4 ppm, 



(4.61) 
(4.62) 
(4.63) 
(4.64) 



where < r p > is the average proton radius connected with the charge distribution. 

For the discussion of the contributions of proton polarizability, which have been obtained on the ba- 
sis of the data of the deep- inelastic scattering of polarized electrons on a nucleon target you see Ref. 17£ ] 



Table VIII. 



Year 



1948 
1948 
1952 
1952 
1955 
1955 
1956 
1960 
1960 
1960 
1962 
1962 
1962 
1963 
1963 
1964 
1965 
1965 
1966 
1966 
1966 
1970 



Reference 



[178 



[178 



[179 



[179 



[180 



[180 



[181 



[182 



[182 



[182 



[183 



[184 



[184 



[185 



[186 



[187 



[188 



[189 



[191 



[192 



[193 



Isotope 



H 
D 
H 
D 
H 
D 
H 
H 
D 
T 
H 
H 
T 
H 
H 
H 
H 
H 
H 
H 
H 
H 



AE, kHz 



1420 410(6) 

327384(3) 
1420 405.1(2) 

327384.24(8) 
1420 405.73(5) 

327384.302(30) 
1420 405.80(6) 
1420405.726(30) 

327384.349(5) 
1516 701.396(30) 
1420 405.762(4) 
1420 405.7491(60) 
1516 701.4768(60) 
1420 405.751827(20) 
1420 405.751800(28) 
1420 405.751827(20) 
1420 405.751778(16) 
1420 405.751785(16) 
1420405.751781(16) 
1420 405.7517860(46) 
1420 405.7517864(17) 
1420 405.751766 7(9) 



Error, ppb 



4224 
9164 
141 
244 
35 
92 
42 
21 
15 
20 
2.8 
4.2 
4.0 
0.014 
0.019 
0.014 
0.011 
0.011 
0.011 
0.003 
0.001 
0.0006 



29 New experiment at the LANL, Los Alamos, is expected to improve the experimental precision by a 
factor 5, cited in [p8[. 



30 It is necessary, of course, to make the corresponding substitutions, 
31 The following results are presented in relative units of Ep. 
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and 
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The published experimental values of the-jpeasurings of the HFS of the ground state in a hydrogen 
and a deuterium are presented in Table VIIll23. 

The difference of values between theory and experiment can be written as 



AE the 



AE exp 



E F 



= (—0.48 ± 0.56 ± uncalculated terms) ppm. 



(4.66) 



The uncertainty 0.56 ppm arises because of the error in the experimental value of the fine structure 
constant and, mostly, the inaccuracy of the data of the proton elastic formfactor. The corrections uncal- 



culated so far could contribute about lppm, Refs. [166, 194 



4.2.4 Muonic helium atom 

The muonic helium atom, which is 4 _ffe ++ ^~e~, has been experimentally discovered in 1975, 
Refs. [195, 196]. From the point of view at a number of electrons this system can be considered as 

which has a middle size be- 
The first measurings of the 



a heavy isotope of a hydrogen having the " pseudo- nucleus" ( 4 He ++ /i ) + 
tween the intrinsic nuclear size and the intrinsic atomic one (~ 130 fm) 
HFS, carried out in 1980, Refs. jl97|-||l99t[, led to the following results: 



AE e ff^z = 4 464.95(6) MHz (13ppm), 



AE e 



and 



AE' 



[199 



4464.02(10) MHz (22 ppm), 



4 465.004(29) MHz (6.5 ppm). 



(4.67) 
(4.68) 

(4.69) 



The last experimental result gave the opportunity to find the magnetic moment of a negative charged 
i nuon 

3.183 28(15) (47 ppm), (4.70) 



Hp 



what allowed one to check the predictions of the CPT- invariance. According to it, the magnetic moments 
of a particle and an anti-particle are to be equal. This quantity for a positive charged muon, which was 



measured more accurately in the experiments of the muonium HFS, Ref. [150 1 , 



^ = 3.183 3461(11) (0.36ppm), 
Pp 



(4.71) 

and, also, in the experiments on observations of the muon spin rotation in liquid, the Larmour precession, 

3.183 3441(17) (0.53 ppm), 



Ref. P00|| : 



Hp 



(4.72) 



is in agreement with (4.7C) to the precision of several tens of ppm. For comparison, the agreement 
between the electron and positron magnetic moments is 0.13 ppm, Ref. [ 201 1 , and is 7500 ppm between 
the proton and anti- proton magnetic moments, Ref. [ 202 1 . 

Theoretical description of a muonic helium atom, which is generally connected with the theoretical 
methods used for muonium, can be found in Refs. [202 [-[206]. We should like to note essential contribu- 
tions of relativistic and radiative corrections to the Fermi energy in this system (M is the " pseudo- nucleus" 
mass), 



E F = l^ a 2 i?ooC !^(i + = 4 5i 6 .9 6 MHz. 

3 m/ M' 



(4.73) 



:S2 



The only experimental result of the HFS of the n = 2 state ( n is the principal quantum number) in 



a hydrogen, known to us, is AE\ 
Ref. HI, 



hfs 



(2S, H) 



177 556.6(3) kHz, Ref. 177 . It satisfies the Breit formula, 
AE(1S) 



AE(2S) = 



5a 2 )' 



(4.65) 
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The obtained theoretical values are I 



AE M$l = 4465.1(1.0) MHz, (4.74) 



, 203[ 

A£*te = 4462.6(3.0) MHz, (4.75) 

A£rf |o| = 4 464.8(5) Mi/z, (4.76) 

A£*te| = 4460 MHz. (4.77) 



, 206[ 

The first value has been obtained by means of the variational methods; the second and third ones, on 
the basis of direct calculation of the Feynman diagrams by perturbation theory method; the fourth, in 
the framework of the Born-Oppenhcimer theory. 



4.3 Fine structure 
4.3.1 Positronium and muonium 

For the first time the fine structure interval has been investigated in a positronium in []20£ 

A£(2 3 5i -2 3 P 2 ) = 8628.4 ± 2.8 MH z (1000 ppm). (4.78) 

In contrast with a hydrogen atom, the mentioned levels are not degenerated in the order of a 2 Roa. 

The experiments of this kind permit checking the validity of QED for the excited states of pure 
leptonic systems. 

The recent experiments [210 achieved the accuracy of 300ppmo 



A£(2 3 Si -2 3 P 2 ) 
AE(2 3 S 1 -2 3 P 1 ] 
AE(2 3 S 1 - 2 3 P ) 



8 619.6(2.7)(0.9) MHz, 
13001.3(3.9)(0.9) MHz, 
18504.1(10.0)(1.7) MHz. 



The theoretical predictions for the first excited states, based on the BS equation, are, Refs. 
IT6| , |i"T7| , [2lH (see also §l|, EDJ)> 

3 



£(l d Si) 



E(2 3 S 1 ) 



E{2 3 P 2 ) 



E{2 3 P 1 ) 



E{2 3 P ) = 



R 

4 



1 

2 

logR(l,0) 



49 
96' 



2n 



a log a 



a 
n 



+ Aisaloga 1 + E>iga 



:R n 



-logR(2,0) 



65 9 3 id 

or H a 3 loqa i H 

192 2tt y 7T 

+ A 2 soi A log a -1 



1 4, „ 

h -loq2- 

15 3 y 



3 r 



97 1, „ 

1- -loq2- 

120 6 y 



B 2 sa 



Ir 



-1 



43 



2 a 

960° ~ ~ 
A P2 a log a -1 + B P2 + . . 



1 



-1 



47 



a" 

or 

192 7T 

Apia^log a -1 + -Bpia 4 
! D / 1 95 2 a a 



' 1 

45 

}. 

_5_ 

36 



-logR{2,\) 



-logR{2,\) 



+ 



Apoa log a 



Bpoa 



+ ...}, 

25 4 
72 + 3 

■•}. 



logR(2,l) 



(4.79) 
(4.80) 
(4.81) 



(4.82) 



(4.83) 



(4.84) 



(4.85) 



(4.86) 



33 r 



3 The calculation of some of the contributions can also be found in [207, 208], but the numerical results 
shown there are highly different from the experimental results. 

34 The first error in these expressions is statistical, the second one is systematical. 
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where logR(n,l) is the Bethe logarithm [214, |215| : 

logR{l,0) ~ 2.984128 5, 

logR(2,0) ~ 2.811769 9, 

logR(2,l) ~ -0.030 016 7. 



(4.87) 
(4.88) 
(4.89) 



The coefficientSpilis and A25 have been calculated recently, but the formulae obtained in |47], Q differ 
from each otherEil 



AE 
AE 



}(a 6 loga) 
} {a 6 loga) 



5 67 -l^/O^sl 

—ma log a 3—. 

24 n J 

1 67 -l^O^sl 

—ma log a 5 — . 

12 J n 3 



(4.90) 
(4.91) 



The coefficients B are not yet calculated. The coefficient 1 of the term aHog a 1 i?/8 gives the contribution 
5.7 MHz and of the term a 4 R/8 gives 1.2 MHz. 

There is a simple formula for the a 2 Roo contributions to the S- levels of an electron-positron sys- 



tem [212 



E(a 2 R 00 ) = 



11 1 

64n 



7_ 

12 



Sit 



(4.92) 



The contributions of the order a 3 R and higher arise from the allowance for the radiative corrections, a 
vertex function, a vacuum polarization, an electron-positron self-energy and an annihilation interaction 
channel. Probably, the formula for the S- states of the electron-positron system, which has been derived 



in Ref. [ 2 1 2 1 , is not correct, as is shown, e. g., in Refs. Eq, 213]. It can be seen from the result of the last 
paper that these corrections are not proportional to 1/nr: 



A£(a 3 i?oo) 




14, 

—log 2 — 6 log a 

O 



3 -logR(n,0)-4(—+log2)8i S 



(4.93) 



The numerical values of theoretical predictions, with taking into account the logarithmic corrections, 
are the following, Refs. |47|, |48|E3: 



AE(2 3 Si 



AE(2 3 Si 
AE{2 3 S 1 

>3< 



1233 607 211.7 MHz 
8 627.7 MHz 
13 013.3 MHz 
18 498.5 MHz 



1 233 6Q7 221.69 MHz, 
8626.21MHz, 
13011.86MHz, 
18497.10 MHz. 



l d Si) 
2 3 P 2 ) 
2 3 Pi) 

AE(2 3 S 1 - 2 3 P ) 

Due to the development of the experimental methods based on the Doppler-free two-photon spec- 
troscopy, now it is possible to measure the "gross structure" interval (IS — 2S) at the accuracy level of 
some MHz. These experiments give the opportunity to find the value of the fundamental constant as 
the Rydberg constant most precisely, see below. 

In Ref. [216 1, the latest results of the measurings of the AE(1 3 Si — 2 3 Si) interval are presented for 
a positronium and a muoniumE3 



AEpf(2 3 Si -l d Si) = 1233 607 218.9 ± 10.7 MHz = -cR^ -83 516.6 ±10.7 MHz, 



(4.98) 



35 To check the results of Fell and Khriplovich et al. the experiments at the accuracy level of ~ 10 ppm 

and the second one, the result of Fell. 



are necessary. 

36 The first column contains the result of Khriplovich et 

37 The results of previous experiments are 



AE 



exjL. 



217| 



41.4(5) GHz; 



(4.94) 
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AEl X jP(2 3 Si - l 3 5i) = 2 455 527 936 ± 120 ± 140 MHz. (4.99) 
In the case of a muonium, there is an agreement with theory, Ref . |222| , 

AE < j^ e u or (2 3 S 1 - l 3 5i) = 2 455 527 959.6(3.6) MHz. (4.100) 

In the case of a positronium the value is 16.6(10.7) MHz greater than the Fulton's theoretical result. 



However, since many authors consider the formula for 5- states in [ 212 1 to be incorrect, it is preferable 
to compare the experimental value (4.9S) with the value of Ref. [213], where the fine structure has been 
found by using a simple potential method, Ref. [223], 



AE theor 

AE^ e u or {2 A S 



{2 3 S 1 ~l 3 S l ) = 
3q - 1 3 Si) = 



83 5Q7AMHz, 
2 455 528 055 MHz, 



(4.101) 
(4.102) 



or with the latest results of Fell or Khriplovich et al. 

The fine structure of the first excited state (n = 2) in a muonium has been investigated in the papers 



[22J-[g26[, see below Subsection [4.5.2 
Let us mention the works 



227], in which the problem of validity of the CP- invariance in the 



lepton sector has been under consideration. The experimental and theoretical limits were obtained there 
for transitions like 2 3 5i — > 2 1 P%. 



4.3.2 Hydrogen and Deuterium 

The results of measurings of the 15 — 25 interval have been reported in [228|-[236| 3 j.The impor- 
tance of .these investigations is clear, because they give information about the Lamb shift value of the 
15- leveltj, which is impossible to find by means of the radiofrequency spectroscopy method used in the 
experiments to find the n = 2 Lamb shift. They also gave a possibility to determine the numerical value 
of the Rydberg constant with the highest precision. 



In Table IX the results of the measurings of the 15 - 
The energy characteristics are presented in MHz. 



25 interval and 15- level Lamb shift are given. 



AE l 
AE f 



21fi| 



1 219| 



1 233 607 185(15) MHz = -cR c 

8 

1233 607142.9(10.7) MHz, 



83 545(15) MHz; 



(4.95) 
(4.96) 



for a positronium; and 



4 AA [22C| 



613 881 924 ± 30 ± 35 MHz, 



(4.97) 



for a muonium (sec for the discussion Ref. 1 22 1 



38 See Refs. [237, 238 for the review of early results of investigation of the interval 2 2 S\/2 — 2 2 P- 



3/2- 



39 The first successful attempt had been made, Ref. [ 239 1 , as early as 1955 to measure this value 
E exp (LS, 15) = 7.9± 1.1 GHz. See also [240], where the value of the isotopic splitting has been discussed 



in a hydrogen for the first time as a result of observations the Layman lines in a hydrogen and in a 
deuterium. 
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Table 1X0 



Year I 


Lefcrencc 


3 Isotope 


AE(1S-2S) 


E(LS, IS) 


1975 


|228| 


D 


- 


8 300(300) 


1975 


[228| 


H 


- 


8 600(800) 


1975 


[229| 


D 


_ 


8 250(110) 


1975 


[229| 


H 


- 


8 200(100) 


1980 


[230| 


D 




8 177(30) 


1980 


[230| 


H 




8151(30) 


1986 


[231| 


H 


2 466 061395.6(4.8) 


8184.8(5.4) 


1986 


[232| 


H 


2 466 061397(25) 


8182(25) 


1987 


[233| 


H 


2466 061413.8(1.5) 


8173.3(1.7) 


1989 


[234| 


H 


2466 061413.19(1.75) 


8173.9(1.9) 


1989 


[235| 


D 


2 466 732 408.5(7) 


8183.7(6) 


1989 


[235| 


H 


2 466 061414.1(8) 


8172.6(7) 


1990 


[236| 


H 


2 466 061413.182(45) 


8172.804(83) 



The result of Ref. p35| of the IS Lamb shift value has been obtained when the following value of the 
Rydberg constant was used: 

Rao = 109 737.315 714(19) cm -1 . (4.103) 

that is the average value of Refs. 258 25E]. There are some values of the isotopic splitting AE(H — D), 
received on the basis of this technique 



AE e 



AE e mm,{H 



3t 



AE 
AE e & 



D) 
D) 
D) 



670 993(56) MHz, 
670 992.3(6.3) MHz, 
670 994.33(64) M Hz, 
670 994.337(22) MHz. 



These values are to be compared with the theoretical values: 

AE t t^(H-D) = 670 994.39(12) MHz, 

y [241 



AE$ffl(H-D) = 670 994.414(22) MHz 



(4.104) 
(4.105) 
(4.106) 
(4.107) 

(4.108) 
(4.109) 



which are in agreement with the last experimental results within an error. 

The numerical values of the IS- level Lamb shift, obtained through the theoretical calculations, 

" l35| 

8 172.89(9) MHz, 
8 173.03(9) MHz, 
8 184.08(12) MHz. 



Refs. p42]-p48|, you can find in Refs. [ p34 

E t ^^(H,LS,lS) 
E t ^^(H, LS, IS) 
E t t^g(D,LS 1 lS) 



(4.110) 
(4.111) 
(4.112) 
(4.113) 

The proton charge radius is supposed in Ref. [ p35| to be equal to 0.862(12) /to, see Ref. |50|li], and the 
deuteron charge radius used in Ref. [241| is 1.962 7(38) /to, Ref. [251]. 

Provided that the Lamb shift value is known from theory, information about the Rydberg constant 
could be received by means of comparison of the calculated value of the IS" — 2S interval with the 
above-mentioned experimental one, see the next subsection. 



40 The result [ 231 1 shown in the Table has been obtained by using Roo of |255]. Better agreement with 
the theoretical result is achieved when R^ of ]253| ] is used, E(LS, IS) = 8 174.8(8.7) MHz. 

41 The previous measuring of the proton radius, Ref. [ 249 1 , should not be ignored because it allows one 
to get a better agreement with theoretical predictions for some experiments. 
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4.4 Rydberg constant 

The latest measurings of the Rydberg constant are the following 
Table X. 



Year 



Reference 



R, cm 



Interval 



1974 
1978 
1980 
1981 
1986 
1986 
1986 
1986 
1987 
1987 
1989 
1989 
1989 
1992 
1992 



252] 



253] 



254] 



231] 



232] 



256] 



257] 



258] 



233] 



259] 



234] 



235] 



236] 



260] 



109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 
109 737 



3143(10) 

314 76(32) 
31513(85) 
31521(11) 
31492(22) 
3150(11) 

315 69(7) 
315 69(6) 
315 73(3) 
315 71(7) 
315 709(18) 
315 69(8) 
315 73(3) 
315 6841(42) 
315 683 0(31) 



2P-3D 
2S-3P 
2S-3P, 2P-3D 
2S-3P 
1S-2S 
1S-2S 
2S-3P 
2S-8D, 10D 
2S-4P 
1S-2S 
2S-8D,10D 
1S-2S 
1S-2S 
1S-2S 
2S-8S, 8D 



12D 



4.5 Lamb shift 

4.5.1 Hydrogen 



In Ref. [ p61[ the results of optical measurings of the Lamb shift of the IS*- level in a hydrogen have 
been reported. 

E(LS, 15, H) = 8 172.82(11) MHz (Uppm). (4.114) 

The technique is based on comparison of frequencies of the two-photon transitions between the IS" — 25 
and 25 — 45, AD levels. It highly differs from the experimental technique of an indirect measuring of this 
quantity, see the preceding subsection, by the two-photon Doppler-free spectroscopy methods in 15 — 25 
transitions. The fi rst optical measuring of the Lamb shift of the 45 level in a hydrogen has also been 
given in Ref. |26l| : 

E(LS,iS,H) = 131.66(4) MHz (300 ppm). (4.115) 
This value could be compared with theoretical predictions of |238[] : 

£||ff|(£5,45, H) = 133.084(1) MHz, (4.116) 
E$$$(LS,4S,D) = 133.254(3) MHz, (4.117) 



and with the radio-frequency measurings, Refs. [ 262 ]-[ 264 1 . 

E e ^^(LS,AS,D) = 133(10) MHz, 
e: ^(LS,4S,H) = 133.18(59) MHz, 
(LS,4S,H) = 132.53t a f 8 MHz. 



E 



E 



1 264 



(4.118) 
(4.119) 
(4.120) 



The discussion concerning the n = 3 Lamb shift in a hydrogen atom can be found in Ref. [265|. Below 
we reproduce the results presented in this article: 



E e: { ^(LS,3S,D) = 314.93(40) MHz, 



(4.121) 
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35, if) 


= 313.6(2.9) MHz, 


(4.122) 


a p66 


(L5, 35, D) 


= 315.3(8) MHz, 


(4.123) 


E §67| 


(L5, 35, if) 


= 315.11(89) MHz, 


(4.124) 




(L5, 35, if) 


= 314.819(48) MHz, 


(4.125) 




(LS, 35, if) 


= 314.898(3) MHz. 


(4.126) 



The precision of these experiments approaches the precision of measurings of the n = 2 Lamb shift, 
which were of great importance for checking the predictions of QED. In Table XI the results of all the 
experiments for n — 2 level are given. 



Table XI. 


Year 


Reference 


AE, MHz 


Error, ppm 


1953 
1969 
1970 
1975 
1979 
1981 
1982 
1983 


I 
I 
I 
I 
I 
I 
I 
I 


268 
269 
270 
271 
272 
273 
274 
275 




1057.774(100) 

1057.772(63) 

1057.90(6) 

1057.892(20) 

1057.862(20) 

1057.845(9) 

1057.8594(19) 

1057.851(2) 


94.5 
59.6 
56.7 
18.9 
18.9 
8.5 
1.8 
1.9 



The above results depend essentially on the parameter r, the life time of the 2P stat e . Si nce the 
experimental data of this constant are absent, the value was found theoretically, Refs. [274, 275 . After 
taking into account relativistic corrections one has 



. /2 xR „ 2 , 9, / 2 xR me 4 a 3 (l + a 2 Zog|) 
7rel . = ATTc(-fR H a 3 (l + a 2 log-) = (-f— 3 i ; - y * J , 



(4.127) 



(Rh is the Rydberg constant with allowance for the finite mass of a proton.). 

The leading radiative corrections, the sclf-cncrgy and the vacuum polarization, contribute additionally 



Trod. 



^c(-) 8 R H a 3 



#(2,1) , s 1 
8 or 



1 

64 



19 
30 



(4.128) 



R(n, I) is the Bethe logarithm. 
Then, we have the numbers 



1 

7 = - 

T 



6.264 8812(20) x 10 s c 
1.596199 46(48) x 10" 



8 „-l 



(4.129) 
(4.130) 



(with taking into account the above-shown corrections). 

In the case of a deuterium (n — 2) we know the following experimental results: 



E e &77^{LS,D) = 1059.00(6) MHz, 



E' 



1059.24(3) MHz. 



(4.131) 
(4.132) 



Let us note that further improvement of the accuracy of the n — 2 Lamb shift experimental values in 
a hydrogen faces the serious difficulties; namely, the natural width of the 2P state is about ~ 100 MHz. 
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The total theoretical formula for the Lamb shift (n = 2) in a hydrogen is, Ref. [277 

1 1 • " " { 1 m 

8 7* 



AE LS = AE- 



2Si 



/2 



AE- 



2 Pi 



a(Za) m , u . » 
- 3 

D7T TO 



log(Za)- 2 - 2.207 909- 



427 3 
nZa(—--log2) + (Za) 2 



~log 2 (Za)- 2 



0,323, 



55 

(4Zo 5 2 + -); 05 (Z«)- 2 



+ {Za) 2 \G s . e .{Za) + G„. p .(2a)] + a(-2 ) ^ + 



(Za) 5 TO 2 
6ttM 



jlog(Za) 2 



12 K ' p 48 

(-0.415 ±0,004)], 



2.399 77 

1 (Za) 4 m 3 



3 7 

-TrZa 

4 



- + Zog(2Za)- 1 -4,25 



M 2 



a(Z«) 5 ra 2 
8M 



,35, „ 39 31 . 

(—Zoo 2 1 )- 

V 4 y 5 192 ; 



(4.133) 



where the self-energy and vac uum polarization contributions (G s . e . and G v . p .) can be expanded in the 
Wichmann-Kroll form §4§]-c,|27|]: 



Gv.p. 



2lS + ^( Za ) l °9iZa)- 2 + 0.5(Za) 



(4.134) 
(4.135) 



G s . e . = -24.1 + 7.5(Za)log(Za)- 2 + 12.3(Za) ± 1.2. 

They give —24.0 ± 1.2, Ref. [ p77| , as a sum in the case of a hydrogen atomElS. 

The numerical value of the Lamb shift of the n = 2 le vel in a hydrogen, corrected with taking into 
account the new calculated corrections, was given in [277 



E,$fSt(LS,H) = 1057.855 ± 0.011, MHz when < r p >= 0.805(11)), 
E$fS(LS,H) = 1057.873 ± 0.011 MHz, when < r p >= 0.862(12)). 



(4.136) 
(4.137) 



For the earlier theoretical works see Refs. [117]-c,[174],[282]-[284], where the recoil corrections of the 

M 



order (Za) 5 ^ have been calculated. The corrections obtained from the diagrams of radiative exchanges, 



which have the order ~ a(Za) 4 ^-, ~ a(Za) 5 ^-, have been calculated in the external field approximation 



in Ref. ]247| |. The corrections of the order ~ (Za) 4 -^- can be found in Refs. 226, 283 1. The contributions 
arising after taking into account the finite size of a proton, have been discussed in J277| , where the 
corrections of the order ~ (Za) 6 ^- have been found out. 

The correction of the order a 2 (Za) 5 m, which. is the binding correction from the two-loop radiative 
exchange diagrams, has not yet been calculatedE 4 ! It is extremely desirable to calculate it in order to 
achieve the theoretical precision ~ 1 kHz. 



4.5.2 Muonium 

At present, two experimental results are known for the Lamb shift 2S , 1 / 2 ~2P 1 / 2 j J — 1 m a muonium 



E\ 



f^(LS,Mu) 



1 0701^5 ± 2 MHz, 



(4.138) 



42 For the case of other ions (Z ^ 1) see Refs. j^g p79 |. 

43 On the basis of the new analy tical method relied on division into the low and high energy part, 
K. Pachucki calculated, Ref. J280[ , the principal contribution to G(s.e.). For the coefficient Ago he 
presented Aqo(1S) = —30.92890(1) and Aqq(2S) — —31.84047(1), which are much more accurate than 
the previous calculations, Ref. J243|, |244j. |246|, |28l|, and the approximation of P. J. Mohr, Ref. |242|-b. 



44 See Ref. [ 285 1 , devoted to calculation of the corrections of this order from the diagrams with polar- 
ization insertions into external Coulomb legs and from the diagrams with the radiative insertions into an 
electron line and one polarization insertion into a Coulomb leg. The calculation of the remained diagrams 
is in progress by this group. 
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lfh(LS, Mu) 



1 054 ±22 MHz. 



(4.139) 



The theoretical value was given by Owen, Refs. [ 224 ]- [ 22G | : 

Ete(LS,Mu) = 1 047.03 MHz, 



as well as the fine structure interval: 2P- 



3/2 



2P 1 



/2 



ffffl(F5,Mu) = 10 921.50 MHz. 



(4.140) 



(4.141) 



In contrast with a hydrogen atom, we have not come across with the problems of structure in a 
muonium, analogously to the calculations of the HFS in such a system. 



4.5.3 Helium 



The precision of experimental measuring, Refs. [288, 289 1, of the frequencies of transitions between 
Rydberg states of A He + reached such values by the beam-foil spectroscopy methods, which made it 
possible to check the Lamb shift value (2Si/ 2 — 2P1/2) m a helium within the corrections of the order 
~ a(Za) 6 mc 2 . In Table XII all the experimental results are presented concerning t he m easurings of this 
quantity, including the early ones based on the Lamb - Rutherford technique, Ref . ||268| . 

Table XII. 



Year 



1950 
1952 
1955 
1957 
1971 
1979 
1987 
1988 



Reference 



[290 1 



[291 



[292 1 



[293 1 



[294 1 



[288 1 



[289 1 



E(LS), MHz 



14 020(100) 

14 021(60) 

14043(13) 

14040.2(1.8) 

14046.2(1.2) 

14040.9(2.9) 

14041.9(1.5) 

14042.22(35) 



Error, ppm 



7130 

4280 
930 
128 

85 
207 
107 

25 



The first theoretical investigations of the Lamb shift in a helium atom have been carried outj-in 
Ref. §97|. The modern calculations, Refs. [§37|, ||| |7|, p89| , give the different results (see also f24§) 



E'^jLS, A He+) 
The latest calculations of G s . e .(Z — 2) give 



G% 



\(Z = 2) = 
i(Z = 2) = 



14 044.5(5.2) MHz, 
14 045.12(55) MHz, 
14 042.36(55) MHz, 
14 042.26(50) MHz. 



-22.8 ±2.0, 
-22.0 ±0.3. 



(4.142) 
(4.143) 
(4.144) 
(4.145) 



(4.146) 
(4.147) 



Besides, we would like to menti on the recent works devoted to investigation of highly excited states of 
a hel ium atom, Refs. [|299| -| 304 . The Table of the latest results for i ons of other atoms was shown in 
Ref. [ 289 1 . The investigation in a muonic helium can be found in Ref. fl305| . 



45 The calculations of Ref. [28£] have been fulfilled by using the new value of the nucleus radius 
1, 673(1) fm, Ref. [|296||, and the obtained result is in excellent agreement with experiment. 
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4.6 The anomalous magnetic moment (AMM) 
4.6.1 Electron 

The opportunity of calculations of the AMM of an electron is guaranteed by the renormalizability 
of QED by means of the expansion in the perturbation series in 2- with the finite coefficients a^. 



o — 2 ot cy. cy 

— = a e (QED) = a u - + a IV (-f + a VI (-f + a V iii(-f + ■■■ 

At present, the value of AMM of an electron is known up to the eighth order. In thc_papers |30 
the calculations of this order have been finished. The obtained result is the followingcB: 



(4.148) 
1C 



theor 
a e,VIII 



■1.434(138). 



(4.149) 



If we use the most precise value for a, the fine structure constant, Ref. [ 3 12 1 , defined by means of the 
quantum Hall effects, 

a' 1 = 137.035 9979(32) (0.02Appm), (4.150) 
then the most precise value for the AMM of an electron i E 

a theor = i 159652 140(27.1)(5.3)(4.1)x 10" 12 



that is in agreements with the 
1.7 standard deviation, Ref. [ 149 1 



(4.151) 

erimental values for an electron and a positron with the precision of 



a ex P = ! 159 6 52 188.4(4.3) x 10" 12 , (4.152) 
a e x + p = 1 159 652187.9(4.3) x 10" 12 . (4.153) 
The theoretical result includes the analytically calculated contributions of the 2nd and 4th orders: 

ajp = 0.5, (4.154) 

_2 



a $3 - 



197 
144 

1 . m 



7i~ 

i~2 



n 2 log2 3C(3) 



= -0.328 478 965 (footnotM), (4.155) 
7 (footnotM), (4.156) 

Lfl llifj, \ lll^ "I'll,/ 

and the contribution of the sixth order, Refs. ]322fl - ||325[ , consisting of 72 diagrams (with not all of them 
being calculated analytically), which have been corrected in Ref. []31( 



[321|,m e , 1 ,m e 2 / m e 4 to 

a e t 7V J ( ) = ) + ( ) l °9 = 5 - 198 X 10 

TO,,, 45 TO,,, \ TO,,, TO 



theor 
a e, VI 



1 



.17611(42) (footnotM), 



(4.157) 



4(i 



The presented result is much more accurate than the preliminary one, = (—0.8 ± 2.5). 

47 The most considerable uncertainty (27.1) arises from the uncertainty of the fine structure constant; 
the second one, from the numerical calculations of a* w ; the third one, from the numerical calculations 

n f theor 
01 "e, VIII- 

48 The hypothesis of Ref. [ |313|1 , that this disagreement (l-7tr) is caused by existence of a scalar electron 
(supersymmetric particle), but not only by the experimental error of a and the numerical integration 
inaccuracy, is not of the present interest. 

49 The history of experimental results is presented in [|314| ]. 

50 Recently the contribution a e jv was recalculated [317 in the Fried- Yennie gauge. The result agrees 
with the previous ones [ 316 1 , [ 318 1-[ 320 1 . 

51 The term of the fourth order a e ,jv(^ £ -), caused by existence of the diagrams with the radiative 
insertion of r- lepton into a vertex, is (-^f-) 2 times smaller ( 4.156 ). The term a &: vi(\ 
(^) 3 (^) 2 (^p) 2 is neglectible small for the modern level of experimental accuracy. 

52 Thc Samuel's result 1 326 , a e> vi = 1> 184(5), is now considered to be overestimated. 
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as well as the small corrected terms appearing as a result of taking into account the muon, r- lepton and 
hadron vacuum polarization loops and the contribution of a weak interaction 



Aa e {muon) = 2.804 • 1(T 12 , 

Aa e (t - lepton) = 0.010 • lCT 12 , 

Aa e (hadron) = 1.6(2) ■ 10~ 12 , 

Aa e (weak int.) = 0.05 • 10 -12 . 



(4.158) 
(4.159) 
(4.160) 
(4.161) 



The recent analytically calculated contribution of the light-to-light subdiagrams in the sixth order of 
the AMM, a e , is of the present interest, Ref. [327]: 



theor 
a e, VI 



ill) 



°4 



6 



C(5) - -tt 2 C(3) 



41 r 
540 ? 



-yiog 2 2 + -log 4 2 



+ 16a 4 - -C(3) - 24ir 2 log 2 ■ 
3 



931 
"54 ' 



0.371005 2921 



oo ^ 

Y — 

i. I On, 



2 n n 4 



0.517 479 061 



what agrees with the numerical estimations 

|3lT 



V VI 



ill) 
(77) 
(77) 



0.36(4), 
0.371 12(8), 
0.370 986(20), 



but disagrees witl 



a e W? (77) = 0.398(5). 



(4.162) 
(4.163) 

(4.164) 

(4.165) 
(4.166) 
(4.167) 

(4.168) 



4.6.2 Muon 

The AMM of a lepton of mass nri\ can be written in the following most general form: 

ai = ai + a 2 ( ) + a 2 ( ) + a 3 ( , ), 

TO2 7713 n%2 rn% 

where m,2 and m 3 are masses of other leptons. 
Like for the AMM of an electron one has 

ai = ai,ii(—) + a l j V (-) + a hVI ( — ) +a it viii(—) +■■■ 



(4.169) 



(4.170) 



It is clear that a2,n — (13,11 — 03, iv — because the Feynman diagrams, which could contribute to 
these terms, are absent. 

Let us consider the modern status of investigations of the muon AMM. The value a\ is known from 
the calculations of the electron AMM. The value a2,iv(-^f-) is not so small in contr ast w ith the electron 
AMM because of the large quantity m^/m e . The result is known analytically, Ref. 33C] 



a2,iv( — ) = ^log—^- — + - 4( ) log— ^ + 3( ) + 0( Y) 

36 4 m e m M m M 



3 m e 
= 1.094 259 6. 



From the above formula the contribution 



2, IV 



( — ) ~ 7.807(5) • 10" 



(4.171) 
(4.172) 



5:5 



Recently, this result has been superseded by the authors, Ref. [ 339 1 . 
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can be deduced, Ref. [331], what gives the contribution to the a M equal to 421.2(3) x 10 2 . 

The term of the sixth order a,2,vi{-^-) which arises from 18 Feynman diagrams, containing vacuum 
polarization loops, is also known analytically, Refs. J321|,[332|-[]334| 



<v/&«-p-)=«rvp-)+A 3 



m e 
3C(3) + 3c 4 



^-| C (2) + 10C(2)/ Off 2- 



31 
27 



k(2) - 4C(2)/ «?2 + C(3) log 



2 2 
9 m e 



ON 



1.944 04(-) d 

7T 



where 



11 j ^2t2 ^ 7 4 ^ 

C4 = — — 7r — — 7r loq 2 — —too 2 — —04. 
648 27 y 27 J 9 

The corrections of the order 0(m e /m^) are known from Refs. [334, 335 ] 



(4.173) 
(4.174) 
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10 , 411, ^ 
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1061 
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ir 2 log2 
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7T H 

108 9 
274511 



TT 2 log 2 2- 



54000 



rrifj, 



(4.175) 



As a result the numerical value is VI 



1.920 4550(2). 

For the first time, the contribution of 6 diagrams containing the light-to-light subdiagrams has nu- 
merically been calculated in Refs. [329, 337, |336| . The analytical result has been obtained recently, 
Rcf. p3§: 
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283' 
12" 



+ (— ) 3 
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9 TO e 



11 



+ (^) 4 

TO M 
13283 



7, ^,to u , 41, 9 / 13 9 517\ , ,rn u , 1 > . . 191 , 

-Zo S 3 (-2 + -Zog 2 (-&) + -7T 2 + — % — + 5 C 3) + — tt 2 

9 TO e 18 TO e \ 9 108 / TO e 2 216 



2592 



0((- 



e ^5^ 



20.9479242(9). 



(4.176) 



Thus, the latest evaluations of this order of the muon AMM give the following numbers, respec- 



tivelyRefs. [329, 335, 338 1 



uaJ 
10 



a% iVI (^;v.p.) = 1.9200(14), 1.92045(5), 1.920455 0(2); 
o£' VJ (f^;77) = 20.9471(29) 20.9469(18), 20.947 924 2(9); 
a^ VI (^;sum) = 22.8671(33), 22.8674(18), 22.868 379 2(11). 

The term of the eighth order to the muon AMM was calculated from 469 Feynman diagrams, each 

54 The result of Ref. ]336| , dyii'll) = 21.32(5), has been invalidated by Samuel in Refs. ]335| , |339| , thus 
finishing the discussion with Kinoshita, see above. 
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containing electron loops of the type of vacuum polarization, or the 
of numerical calculations of the eighth order was presented in 1 34f| 



777 

<vm( — ) = 126.92(41) 

TYlf= 



rto-light subdiagrams. The result 



(4.177) 



In Refs. j 342- 1 344] the asymptotic (when — > 0) contributions to the muon AMM have been 

obtained in the analytical form by using the renormalization group techniqucS. They come from the 
diagrams of the eighth order with one loop of electron vacuum polarization, Ref. Q344], 



^vmM = -v l09 ™ e + 48 + 8 C(2) " C ^ 2 + 

-0.290 987... 

and from the diagrams with two loops of electron vacuum polarization, Ref. [342. 

1 
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^oo,2 
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1 , 2 «V 

— log — — 

12 " 



m e 



;C(3) 



m e 



1531 



1728 12 



^C(2)-^SC(3) = 1-452 570. 



1152 



(4.178) 



(4.179) 



Moreover, the result for the diagram of the tenth order (additional loop of vacuum polarization) is 
now known, Refs. [344, 345 1, 



oo.2 / \ 



+ 



1 , 9 , m Ll \ 

2i l ° 9 ( - ) + 



C(2) 
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3409 
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288 

-1.3314. 
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C(5) + ^C(3)--C(2)^2+ 



C(5) 
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128 



32 



C(3) 



|C(2)^ 5 2-^C(2) 



3456 

The numerical estimation of the tenth order is, Ref. H34i 



(4.180) 



aa,x(— ) = 570(140). 

TO R 



The terms 



to m m M 

03,V/( — , — ) 



0/3, VIIl( 



5.24(1) x 10~ 4 
0.079(3) 



(4.181) 

(4.182) 
(4.183) 



The preliminary value has been given in [337 



*2, VIII 



140(6). 



56 When we have submitted this review to print we learned about Kinoshita's recalculation 341] of his 
result of 1990, namely, of the contribution of the eight-order vertices containing sixth-order one-electron- 
loop vacuum polarization subdiagrams. Including new result, which is close to the asymptotic analytic 



result of Broadhurst et al. [344], the value vin(~£r) — 127.55(41) supersedes the result ( 4.177 ). Thus, 
the improved value of ( jOgg ) is a^QED) = 1 165 846 984 (17)(28) x 10~ 12 



57 When we have submitted this review to print we learned about the remarkable results of the pa- 
per [346] where the [n + 1)- loop contributions to muonic anomaly have been obtained analytically from 
the 1-loop diagra ms w ith insertion of the n- loop photon propagator containing n — 1 electron loops. 



58 See also Ref. [ 347 , where the same technique has been used for calculation of the Kallan-Simanzik 
(3- function in the eighth order. 
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have also been calculated there. Combining the above results with, the term a\ known from 
Kinoshita obtained the pure quantum-electrodynamics contribution^ 

a^QED) = 1165 846 947(46) (28) x 1(T 12 . 



The recalculated result of Samuel, Ref. [335], is 



a^QED) = 1 165 846 950(28)(27) x 10 



-12 



After adding the hadron contribution, Refs. p48]]-p5l|P: 

a^(hadron) = 7.03(19) x 10~ 8 

and the weak interaction one, Refs. m 

a^weakint.) = 181(1) x 10~ n , 
one gets the following value for the muon AMM, Refs. |335| , |340| , respectively: 

116 591902(77) x 10" 11 , 



a 



u. the or 



[335 



u. the or 

r 40| 



116 591920(191) x 10" 



(4.184) 

(4.185) 

(4.186) 

(4.187) 

(4.188) 
(4.189) 



which are in good agreement with the known experimental values, Ref. [354] 

„ex P ^ 1 165 937(12) x 10" 9 , 



exp 



1 165911(11) x 10" 



(4.190) 



Let us point out that the experiments planned at B NL, Brookhaven, will increase twenty times the 
accuracy of experiment (up to ±0.05 x 10 -8 ), Ref. [355|. 



4.6.3 t - lepton 



The theoretical calculations of the r - lepton AMM have been carried out in Refs. [356, 357 . The 
obtained estimations are the following: 



a T = 11773(3) x 10" 7 . 



(4.191) 



Experimentally, the quantity a T is investigated less completely than theoretically. The direct ex- 
periment is lack. The cause of this bad situation is a very small life-time of r- lepton, Ref. |35S] 



T = 2.957 (32) x 10 c. However, recently the experimental methods to check theoretical predictions 



have been proposed, Refs. J359( |- J362| . Moreover, some constraints, Ref. [ 363 1 , have been obtained from 
the analysis of electroweak experimental data: 



- 8 x 10~ 3 < a T < 1 x 10~ 2 (2a) 



(4.192) 



For the completeness let us reproduce the data on the energy levels of t + t~ atom, Ref. 1 364 ] . It can 
be made in e + e~ annihilation below r pair threshold. It is clear that the energy levels are defined from 
the formula: 



E n = — - 



m T c 2 a 2 23.7 



4n 2 



keV, 



(4.193) 



59 The first error is the estimation of the theoretical uncertainty, mostly from a,2,vi{-£f-), the second 
one is from the uncertainty of a determined in quantum Hall effect experiments. 

The recent result of Dubnicka el al. |551| is more exact than the result presented by Kinoshita, 



60 r 



[351| 



61 



(hadron) = 6.986 ±0.042 ±0.016 x 10 8 . The model and experimental errors were also given there. 



The shown uncertainty is caused by the lack of information about Higgs mass. 
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where m T = 1777.8 ± 0.7 ± 1.7 MeV/c 2 , Ref. fl365| ], has been used. The decay widths are 

a 5 m T c 2 1.8 x 1(T 2 



T(t+t- 
T(t+t- 



2 7 ) 
37) 



2n 3 n 
2(tt 2 - 9)a 6 m T c 2 
9n3 



■eV, 
1.7 x 10- 



eV. 



More complete consideration of r- physics you can find in the excellent reviews, Ref. [ 366 1 . 



(4.194) 
(4.195) 



4.7 Fine structure constant 

The experimental values for the fine structure constant are following, 
a) The fine structure constant defined from a quantum Hall effect is, Ref. [ 3 1 2 1 

a^iQHE) = 137.035 997 9(32) (0.024ppm); 



(4.196) 



b) on the basis of calculations of the eighth order of the AMM of an electron and its comparison with 
the experimental values, Refs. [149, 31C|e3, 



or x {a e ) = 137.035 992 22(51)(63)(48) (0.0069 ppm); 



(4.197) 



c) on the basis of Josephson's effects, Ref. [367], 

a" 1 (JE) = 137.035 977 0(77) (0.056ppm); (4.198) 

d) on the basis of comparison of the experimental and theoretical results of the HFS in a muonium |l37| 

a^in-hfs) = 137.035 992(22) (0.16 ppm). (4.199) 

The results (a) and (b) agree with each other up to the error level 0.05 ppm, thus proving the validity 
of QED to the above accuracy. However, further enhancing of experimental accuracy of measuring of 
a~ 1 {QHE) and on the basis of other methods is desirable. 



A possible or igin of the discrepancy between ( 4.196 ) and ( 4.197 ) is a possible subquark structure of 
an electron, Ref. [368], and-the some unknown type of an interaction mediated by particles more massive 
than W ± and Z°- bosonsEJ. 
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62 The first uncertainty originates from the experimental uncertainty, the second one is from the uncer- 
tainty of the sixth order calculations, the third one is from the uncertainty of the eighth order calculations. 



The present result, Ref. | pl0| , is 3 times more precise than the one in 1 312 and, moreover, than the ones 
obtained by the other methods. 

63 However, the contribution, caused by this interaction, can not be large because of the m/M term. 
The value ^ eor is more sensitively to the case when there is an unknown light particle interacting with 
an electron, Ref. |70||. 
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